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A USERS MANUAL 
FOR COMPUTERIZATION 
OF HYDROGEOLOGIC DATA 
FOR A RIVER ALLUVIUM 

Cy James W. Naney, Bill B. Barnes, and Douglas C. Keiit^ 

ABSTRACT 

A system for rapid storage and retrieval of data from a mLiltilayered 
aiiuifer based on tlie hydraulic properties of the aquifer was developed, 
'file manual con.si.sts of the data-i’etrieval system, control-card setups and 
siiocifications, lists of variables, formats and source listings for data stor- 
age and retrieval, and graphical illustrations demonstrating how the system 
works. Data for the examples were collected in the alluvium of the 
W’ashita River near Chickasha, Okla, These data were processed according* 
to the manual to produce a lithologic ci’oss section, an isopachous map, 
a specific subdatuni map, a distributed test-hole-data map, and a perme- 
ability-distribution map. The system proved capable of storing and re- 
trieving a large volume of layered hydrogeologic data, and it will have 
wide g'eographie application where hydrogeologic studies are undertaken. 
The system may also be adapted to other types of materials because it 
uses numerical coding exclusively to de.scribc the hydrog'eolog'ic data to be 
.stored, retrieved, and manipulated. KEYWORDS: computerization of 
hydrogeologic data, ground water, ground-water storage, ground-water 
transmissihility, hydrogeology, hydrology, lithology, soil permeability. 


INTRODUCTION 

A layered lithologic system created the need 
for specialized computer iirograms to store, 
retrieve, and manipulate hydrogeologic data. 
Existing prograni.s (/, o, 7 )- would store data 
on labeled punclied cards or provide for the 
aljihameric coding of data on punched cards 
and store and retrieve ground-water hydro- 
graphs. This manual, however, has been 
developed for the selection and retrieval of 
ground-water data ba.sed on the hydraulic prop- 
ei'lies of an alluvial aquifer. It documents tech- 
niques described in an earlier paper by Kent et 


I Geologist anti compulor specialist, Agi’ieultural Re- 
search Service, U.S. Department of Agriculture, 
Chickasha, Okla. 73018, and associate profe.s.sor of 
geology, Oklahoma Stale University, Stillwater, Okla. 
74074. 

- Italic nurahers in parentheses refer to items in “Lit- 
erature Cited” preceding appendix A. 


ah (./) ; excerpts from that paper are included 
here for convenience. 

The retrieval system was apiilied to define 
boundaries and internal hydraulic character- 
istics of an aquifer in the Washita River al- 
luvium near Chickasha, Okla. (fig. I). The 
alluvium is characterized by several discon- 
tinuous layers of silty clay, sand, and gravel 
similar to those found in stratified alluvial fans, 
as well as some stratified basin and Coastal 
Plain sediments. Although the approach used 
in this .study was specifically designed for al- 
luvial aquifer systems, as shown by Naney and 
Kent ( 6 ), it can be used for other similar 
aquifers. 

The computer codes presented in this manual 
are written in FORTRAN IV for the IBM S60- 
65 computer. Some of the techniques used for 
processing the data in the examples were de- 
veloped on the IBM 1130, A restriction on the 
word length for coding some of the data was 
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Ficure 1, — Hydiaulic-coefficient ranges in the Washita River alluvium. 


encountered with tlie IB'Sl 1130. The adaption 
of tlie>e jirograms to both computers, however, 
lias expanded the availability of the technique 
fur which this manual is written. 

Three types of data are documented and 
st(ired in this data-handling; system: (1) gen- 
eral test-hole information (test-hole data); 
(2) stratification characteristics (layer data); 
and (3) water-level rGCord^ (water-level 
data). Tlie test-hole data consist of the well 
number, location, elevation of ground surface, 
elevation of tojj of pipe (well casing), the num- 
ber of layers, and the method of drilling. Layer 
data include layer miniber and thickness, meth- 
od of analy^is, lithologic type and color, and 


approximate range of the hydraulic coefficient 
for each layer. The water-level data include ^^'"ell 
number, record type, date, time, and water-level 
elevmtion. 

DATA-RETRIEVAL SYSTEM 

The data-retrieval system (fig. 2) is based 
on descriptions of the lithology in terms of the 
type and color of material and of the hydraulic 
properties associated with each layer. These 
descriptions are stored by test-hole number, 
latitude and longitude, and layer numlier on 
a disk or magnetic tape. General test-hole in- 
formation and water-level records are stored 
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Figure 2* — ^Type and function of data processed in computer storage and retrieval system. 








Figure 3. — FIo^^^ chart for retrieval of hydrogeologic data, 













by well number and location on separate disks 
or magnetic tapes. Specific information related 
to each of these data is numerically coded and 
stored by varia))le number. Combinations of 
tliese stored data are used as criteria in the 
selection of other data or as data that are to 
be shown grapliically. 

The hydraulic characteristics of a multi- 
layered aquifer system can be evaluated by 
options that allow selective retrieval of data. 
Ten options for data selection and comparison 
are specified on control cards. A general flow 
cliart of the retrieval program is shown in 
figure 3. Figure 4 shows how each specific 
option is selected. Control-card layouts for 
retrieval options A ilu’ough J and a description 
of card types A through N appear in appendix 
A. The variable name and columns assigned 
to each code arc shown for each option. A list of 
variable names with descriptions also appears 
in appendix A. In addition to the options de- 
scribed, subroutine packages may be introduced 
to be called by specified control cards. 

Datastorage format — ^The data-storage pro- 
gram reads data from punched cards in the 
formats shown in appendix B and stores them 
on magnetic tape. Twenty-one variables are 
stored for each test hole, including variables 
that describe the entire well, such as the type 
of well casing or the elevation of the top of 
the pipe. These variables are presented in table 
1 . 

Twenty lithologic zones or layers are allowed 
for each test hole. Currently, no more than 14 
distinguishable layers have been stored for 
any one test hole. There are eight variables 
describing the lithologic characteristics of each 
layer, such as layer thickness or the type of 
material. These variables (table 1) are read 
from punched cards according to the format 
shown in appendix B and stored on magnetic 
tape. 

Ground-water levels are stored on magnetic 
tape for use in conjunction with the test-hole 
and lithologic-layer data. The format used on 
inmched cards to store the data is shown in 
appendix B* Storage and retrieval programs 
for test-hole, layer, and water-level data, along 
with comments and execution cards, appear in 
appendix C. The three basic graphical forms 
produced are test-hole cross sections, lithology- 
distribution maps, and isopachoiis maps. A 
separate retrieval program (not shown in fig- 


ure 2) provides data for maps that show the 
distribution and change in level of the ground- 
water surface. 

Coding and storage of lithologic data . — The 
lithologic type and color are described with 
unique combinations of numbers (table 2) , Only 
five numbers representing grain size are used in 
the code to describe the lithologic type because 

Table 1. — Data-processing code for tape stor- 
age, entire ivell 


Variable Term Code form’ 


1 Well number or test-hole 

number XXXX 

2 Watershed iiumher XXX 

3 Location (latitude: degrees, 

minutes, seconds) XX XX XX 

4 Location (longitude: degrees, 

minutes, seconds) XX XX XX 

5 Top of pipe elevation * XXXX.XX 

6 Ground-surface elevation XXXX.XX 

7 Bedrock elevation XXXX.XX 

8 Total depth of test-hole log (ft) XXX 

9 Distance from logged well to 

observation well (ft) XXX XXX 

10 Total depth of observation well 

(ft) XXX 

11 Drilling method X 

Core 1 

Rotary 2 

Split spoon 3 

Auger — 4 

12 Casing type X 

Metal 1 

Plastic 2 

13 Casing length (ft) XXX 

14 Casing diameter (in) XX, XX 

16 Screen XXX, XX 

Screen length (ft) XX, XX 

Screen diameter (in) . • . XXX. XX 

10 Pump type X 

Submersible 1 

Jet 2 

Turbine 3 

Centrifugal 4 

17 Pump use X 

Domestic 1 

Irrigation 2 

Experimental 3 

18 Pump capacity (gal/min) •• • XXXX 

19 Number of stages XX 

20 Pump diameter (in) XX.X 

21 Number of layers XX 


^ An X represents a space provided fox' the variable- 
number code, A number indicates a specific descriptive 
adjective to bo coded in the space provided for that 
valuable. For example, a 4 put in the space design a tod 
by the X after “drilling method” would only select wells 
drilled by an augei% 
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READ FftOH CARO TYPE A 


( f rniniJiitei- storage limitations. Tiie codes are 
hy iitljologiL layer m tlie computer. Two 
ioaiiaiy i-r mie primary and one secondary 
I'l'ii'VitHfi.t de.'.ci’ilie ca* li litliologic layer. Tlie.se 
d.iM .lie ;jfau rediued to a form that can show 
tie- iiydraiilii properties (cocfficient.s) of tlie 
.xijlJU'er, 

The liUmk.gic type and color are based on 
log dc'criptions made m Die field. These log.s 
iiee.l to he interpreted accurately before tJie 
data I an be standardized ami stored as layer 
data. I- or analytical interiiretation of test-hole 


Figure 4. — Flow chart for 

logs, samples are taken from holes cored ad- 
jacent to formerly drilled test sites. Grain size 
of each core sample is visually described and 
analyzed for comparison with descriptions 
based on drill cuttings. The median grain size 
is equivalent to the grain size having the high- 
est percentage of occurrence in the drill-cutting 
sample and is termed “primai-y constitutenl.” 

giciin sizG 1 epreseiiting* r smaller pGrcent<ijg6 
of the sample is labeled "secondary constitvi- 
ent." With these comparisons and terms, quanti- 
tative relationships can be transposed into 
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retrieval control cards, A-N, 


terms for coding the visiuil description of drill 
cuttings. 

Hydraulic properties of the aquifer are de- 
termined by using a relationship between grain 
size and the hydraulic coefficients (fig, 5). 
The hydraulic coefficients (the coefficients of 
permeability and storage) are based on field- 
pump tests of selected alluvial material, lab- 
oratory tests of the core samples, and grain 
size-specific yield relationships ( 2 ). The hy- 
draulic properties are coded by dividing the 
primary and secondary grain sizes into four 


identifiable ranges (table 3) . Undifferentiated 
sand is included in a fifth range. 

A computer program (appendix D) is used 
to assign a range number (1-5) to the layers 
of each well on the basis of the coded lithologic 
description shown in table 4. Hydraulic-coef- 
ficient values representing these ranges can be 
arbitrarily selected from a graphical relation- 
ship similar to the one shown in figure 5. The 
grouping technique makes them applicable to 
all drill-log data, with only the values repre- 
sented by each code being subject to change 
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Tadle 2. — Data-processing code for tape storage^ stratigraphic layers 


Variable Term Code form^ 


1 Layer number XX 

2 Layer thickness XX 

3 Type of material:- column ly num- XXXXX 

her of primary constituents (zero 
if unknown) ; cohimn code of 
primal y constituent; column 3, 
code of adjective for primary 
constituent (zero if primary 
constituent has no adjective) ; 
column /f, numeral 1 if second 
primary constituent or secondary 
constituent has no adjective 
(zero if no secondary constitu- 
ent)* or code of secondary con- 
stituent if it has an adjective, or 
code of second primary constitu- 
ent with adjective; column 5, 
code of second primary constitu- 
ent with no adjective or second- 
ary constituent with no adjective 
(zero if no second primary con- 
stituent or secondary constituent 
exists), or code of adjective for 
secondary constituent, or code of 
adjective foi second primary 
constituent. 


Constituents X 

Clay ^ . . . - . . 0 

Silt 1 

Sand 2 

Gravel 3 

Sandstone 4 

Carhonate ... .... 5 

Shale 6 

Uvaporite . . 7 

Undifferentiated alluvium .* 8 

Undifferentiated bedrock 9 


Variable Term Code form^ 


3 Type of material — Continued;- 

Log-description adjective . X 

Very fine 2 

Fine 3 

Medium 4 

Coarse B 

Lime ^ 

Dolomite T 

Halite ^ 

Fossiliferous 9 

4 Color of material: major color is XX 

the left digit, minor color the 
right digit (a zero in either place 
indicates no color). 

5 Analytical method ■ ■ * ■ • X 

Visual accumulation tube 1 

Sieve 2 

Driller inspection (field) . 3 

Microscopic examination . . 4 

Hydhaulic coefficients J 

6 Transmissibility (T) (gpd/ft) ■ • X 

Values from 0 to 8,000 .... 1 

Values from 8,000 to 30,000 ■ • 2 

Values from 30,000 to 80,000 . 3 

Values from 80,000 to 160,000 4 

7 Permeability (K) (gpd/ftU . X 

Values from 0 to 80 - . ■ 1 

Values from 80 to 300 .... 2 

Values from 300 to 800 - • . 3 

Values from 800 to 1,600 ... 4 

8 Storage coefficient (5) (dimension- 

less ratio) X 

Values from 0.0000 to 0.0001 t 

Values from 0.0001 to 0.0009 . . 2 

Values from 0,0009 to 0.0020 . 3 

Values from 0.0020 to 0.0040 • 4 


^ An X rep resents a .space provided for the variable-number code. A number indicates a specific descriptive ad- 
jective to l)e coded in tlie space provided for that variable. 

- For example, a log desciiption of a primary and a secondary constituent as ^hnedium sand, with some fine sand*’ 
results in a code name of **san<i, medium, with some sand, fine,*’ or a number code of 12423. Again, a log descrip- 
tion of two primary constituents as ^‘medium sand, and fine sand’* results in a code name of ^^sand, medium, and 
sand, fine,” or a innnher code of 22423. 

For example, a sample with a major color of red but no minor color would be coded 20, one with no major 
color but a minor color of red, 02; a sample with a major color of red and a minor color of brown would be coded 
26, one with a major color of brown and a minor color of red, 62. A code of 00 indicates no color at all was re- 
corded. 

* Code numbers 1-4 represent a range of values for T, K, and S, 
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based upon the geological materials being 
investigated. 

Retrieval of lUkolofjic data . — Data can be re- 
trle^’ed in many forms, but the most useful are 
cross sections and maps that define aquifer 
boundaries and the horizontal or vertical dis- 
tribution of the hydraulic properties in the 
aquifer. The schematic diagram of the retrieval 
system (fig. 6) shows the assignment and use 
of the ranges and corresponding estimated 
values of the hydraulic coefficients for presen- 
tation in different graphical forms. Retrieval 
and testing of test-hole, layer, and water-level 
data are essential to provide information neces- 
sary for selection and plotting of specified 
ranges of the hydraulic coefficients (fig. 2). 
Retrieval of test-hole data provides information 
necessary for selection and plotting of tlie 
ranges at a specified depth or subdatum eleva- 
tion. Water-level data are used to select ranges 
characteristic of each laj'er within the satu- 
rated zone. Testing on the range number of the 
hydraulic coefficient permits the selection of 
specified ranges of the hydraulic coefficients. 


Table 3. — Lithologic type and associated range 
of the hydraulic coefficients 


Rangel 

Lithologic description of log- 

1 

Silt; - silt with some very fine sand; veiy 
fine sand with some silt; silt and very 
fine sand. 

2 

Very fine sand; very fine sand with some 
fine sand; fine sand with some very 
fine sand; very fine .sand and fine sand. 

3 

Fine sand; fine sand with some medium 
sand; medium sand with some fine 
sand; fine and medium sand. 

4 

Medium sand; coarse sand and gravel; 
medium sand with some coar.sa sand and 
gravel; coarse sand and gravel with 
some medium sand; medium and coarse 
sand. 

5 

Undifferentiated material in ranges 2, 3, 
and 4, 


^ Zero is used when a range cannot be assigned to a 
layer. 

- Combination of primary and secondary constituents 
that are classified as indiented. 

Silt and clay arc undifferentiated. 


Table 4. — List of codes repi^esenting four (/ 7 'oups of relative permeability 


Clay and silt 
( 1 ) 


Very fine to 
fine 
( 2 ) 


Sand 


Fine to medium 
(3) 


Medium to 
coarse 
(4) 


10000 

20011 

12211 

22211 

10011 

20012 

12311 

22311 

10012 

20013 

12011 

21122 

10013 


11122 

21123 


20017 

11123 

22223 

10017 

20019 

12200 

22322 

10010 

23010 

12300 


10125 

21010 

12223 


10124 

21012 

12322 


10123 

21013 




10122 

12210 

12310 

12410 

12610 

12010 

13010 

11000 

11010 

11012 


12324 

22324 

12500 

22524 

12423 

22423 

12524 

22425 

12300 


12425 

23125 

12400 


13000 

23124 

12413 


13012 

22513 



13125 




13124 




12513 



11013 

11017 

11010 

13011 

IGOOO 
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Figuke 5. — RelatioiiHhip between hydiauhc-coefficient ranges, median grain size, and 

the coefficient of permeability. 


One important aspect of the retrieval system 
is the use of water-level records for selecting 
data in the saturated zone. A subroutine which 
computes a weighted average of all measiu’e- 
ments within a selected time period is used to 
determine the ax’erage elevation of the ground- 
water surface in each well. 

GRAPHICAL PRESENTATIONS 

Examples are given of five basic types of 
graiiliical iireseiitations of data retrieved. One 


cross-section display and four possible map- 
view displays are shown. A flow chart of the 
steps followed in preparing graphical displays 
and the list of variable names, with descrip- 
tions, used to plot selected data are shown in 
appendix E. 

Cross sections of drill-hole logs. — Drill-hole 
data are punched on cards and may be used 
as input data for plotter programs. The con- 
trol-card setup required to retrieve data for the 
plotting of cross sections is shown in figure V. 
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Figure 7. — Control cards for retrieval of data for plotting cross sections. 
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Figure ft. — Control cards for retrieval of data for mapping distributions of selected 


hole data. 


SUGAR 



Figure 10 — Distnliulion of total depth of test holes. 


For plotting’ options availtible, a description of 
each retrieval control card used is found in 
appendix A. 

A plotted cross section of drill holes near 
Alex, Okla., is shown in figure 8. Appendix F 
sliows the data cards used to plot cross sections. 
Tlie variable name, type of plot, and column 
in wliich each variable is coded are shown for 
each card. Tliese data may then be used on tlie 
plotting equipment available to the user. 


M<i}}s of distributed drill-hole data . — The 
total depths of the test holes were plotted as a 
distributive map. The control cards for the 
plotter program are listed in figure 9, and a 
map showing the distribution of the total depth 
of each test-hole log appears in figure 10. The 
map shows the depth to which lithologic data 
are available, since some lest holes were not 
drilled to bedrock. Appendix F contains the 
data cards used to plot selected data maps, with 
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Figure 11. — Contiol cards for retrieval of hydraulic-coefficient ranees at a specified 

subdatinn. 



Figure 12.~Map of hydraulic ranges for a specified siihdatiini elevation. 


the type of plot, variable name, and column 
containing the variable code indicated. Any of 
the 21 test-hole variables may be retrieved for 
plotting as a distributive map. Appendix E 
gives a list of variable names and their descrip- 
tions for all example plots. 

Mafis of data from specified elevations . — 
Because hydrogeologic data concern informa- 
tion that must be gathered at some depth below 
the surface of the ground, a convenient method 
of looking at various levels, or subdatum eleva- 


tions, was developed by the authors. Digital 
computer techniques are used to develop a 
planar surface of lithologic data for a specified 
elevation above mean sea level. Figure 11 shows 
the control cards for the subdatum searching. 
A map of the test holes that have lithology with 
hydraulic coefficients ranging from 1 to 5 at 
the specified subdatum is presented in figure 
12. Appendix F contain.s the data cards needed 
to plot the map shown in figure 12. 
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Figure 13. — ContJ’oJ cavtls for retrieval of weighted average.^ of the hydraulic co- 
efficients of all layers in each hole. 



ilfrtjis of average eoeffmente of penncahility 
avd storaf/e . — When the ground- water flow 
through a specified cross section or the amount 
of ground water stored in a region of a water- 
shed is needed, it is sometime.^ desirable to 
average the coefficients of )>evmeability and 
storage of the layers that the test hole pene- 
trates. Figure 13 shows the control cards for 
averaging these hydraulic pararoelevs. 

Those averaged values of K (permeability) , 
when displayed in map view (as shown in fig- 


ure 14), provide a distribution of test holes 
that are potentially high or low in ground-water 
yield. A distributive map of 5 (storage co- 
efficient) may be produced using the system. 
It is possible to determine suitable well field 
locations and possible artificial recharge areas 
using areal maps (fig. 1) developed from the 
distributed test-hole data. 

The retrieved data cards used to plot the 
map shown in figure 14 are presented in 
appendix F. 


W 


499 500 501 * 

ABC 






I 

A a 


CARO TYPE E 


CARD TYPE D 


513 51 4 

0 P 

CARD TYPE C 


CARD TYPE 8 


CARD TYPE A 


Figure 15. — Control caiclb for retiieval of an n>opachous map of lithology having 

peimealjility range b 1-5. 



Figure 16, — Isopachous map of hydiaulic properties, ranges 1-5. 


Isopachous maps of specific hydraulic prop- 
erties . — An isopachous map of all the lithologies 
\vithin an alluvial system having the same 
hydraulic properties was also plotted. Such a 
map is used by planners when well fields are 
being installed or augmented. The amount of 
ground water available in an area may be de- 
termined from such a map (J). 


A control-card setup for retrieving data for 
an isopachous map of all subsurface materials 
having a hydraulic code 1-5 is presented in 
figure 15. The isopachous map plotted is shown 
in figure 16. 

Appendix F shows the retrieved data cards 
used to plot the map shown in figure 16. 
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APPENDIX A.-CONTROL-CARD SETUPS, CARD SPECIFICATIONS, 
AND LIST OF VARIABLES FOR OPTIONS A-N 


Control-Card Setups 




Figure A-2. — Retrieval option B : type of material. 


17 






501 502 503 

ABC 




CARD TYPE Q 


CARD TYPE F 


515 5 16 

0 P 
■ CARD TYPE C 


CARD TYPE B 


CARD TYPE A 


Figure A--6. — ^Retrieval option F: data for logging a test hole. 






















































































Figure A-13. — Description of card types G, H, and I 
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PiOUBE A-14. — Description of card types J, K, and L. 
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Figure A-16.— Description of card types M and N. 


List of Variables 

Variable 

name Option Format Description of variable 

Card Type A 


NWELL 

A-J 

16 

LVAR 

A 

16 

NUM 

A 

16 

LVARl 

C-J 

16 

NVAR 

A-J 

15 

ICOND 

B 

16 

ICOL 

C 

16 

IDAT 

D 

16 


Number of test holes to be searched for data 
retrieval. 

Number of lithologic-layer variables to be 
retrieved. 

Number of selected values of lithologic-layer 
variables to be retrieved. 

Number of additional lithologic-layer data 
to be retrieved. 

Number of descriptive well data to be re- 
trieved. 

Number of lithologic materials to be re- 
trieved. 

Number of lithologic colors to be retrieved. 

If a 1 is punched in column 40, data are 
retrieved from subdatum planes (see card 
type J) . If column 40 is blank, this 
option is bypassed. 
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Vcii mhlc 

mune 

Option 

Format 

Description of variable 

INC 

E 

Caid Typo 

15 

A— Continued 

An integer punched in columns 41^46 is the 

ILOG 

F 

16 

number of inciements to be searched 
above bedrock (see card type K) . If 
columns 41-45 are blank, this option is 
bypassed, 

Alin column 50 indicates a total test-hole- 

ILVL 

G 

15 

log retrieval. If column 50 is blank, only 
the selected layer variable is retrieved. 

A 1 in column 55 calls water level for 

IGSUB 

H 

15 

defining saturated zone (see caid type 
N ) . If column 66 is blank, this option is 
bypassed. 

An integer in columns 5G-G0 indicates the 

KTESK 

I 

15 

number of feet below ground surface to 
be searched by increments. If columns 
56-60 are blank, this option is bypassed. 
An integer in columns 61-65 indicate.s the 

KTESS 

I 

16 

number of permeability values assigned 
(see card type L). If columns 61-65 are 
blank, this option is bypassed. 

An integer in columns GG-70 indicates the 

KSAT 

I 

15 

number of storage-coefficient values as- 
signed (see card type M), If columns 
66-70 are blank, this option is bypassed. 
If column 7B is a 1, JC and S values are 

IPUNII(l) 

D,E,G 

averaged for the total well log. If column 
75 IS blank, K and S values are averaged 
only for tl^e saturated zone. 

Card Type B 

15 A 1 in column 6 causes retrieved layer data 

IPUNH(2) 

J 

16 

to be punched on cards, 

A 1 in column 10 causes the net thickness 

IPUNH(3) 

B--J 

15 

of selected lithologic layers to be punched 
on cards. 

A 1 in column 16 causes test-hole (well) 

IPUNn(4) 

C,P 

16 

variables to be punched on cards, 

A 1 in column 20 causes lithologic color to 

IPUNH(5) 

B 

IB 

be punched on cards. 

A 1 in column 25 causes lithologic material 

IPUNIUG) 

I 

15 

type to be punched on cards. 

A 1 in column 30 causes the weighted aver- 

Note: 

age of hydraulic coefficients to be 
punched on cards, 

A blank in appropriate columns ^vill bypass punch output for that option. 

IWELL(n) 

A-J 

Card Type C 

1616 This array is used to select test-hole data 

IVAR(w) 

A 

by numbers. As many as 1,000 drill-hole 
numbers may be retrieved from this array 
selectively. 

Note: n indicates array size. 

Card Type D 

816 Data are retrieved from this array for as 


many as 8 tithologic-layer variables. 
These variables describe hydraulic and 
geologic characteristics for discrete 
layers. 
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Format 


Variable 

lUime 

TVALU(h) 

lADVR(n) 

ISTOR()0 

MAJTS{n) 

MACOL(n) 

SUBDA 

SUBTS 

SUBIN 

AINC(h) 

AKVAL(m) 

SVAL(?0 


Option 


A 


C-J 


A*-J 


B 


C 


D 

D 

D 


E 


I 


I 


DescrijHion of variable 


Card Type E 

6P12,4 Data are retrieved from this array for as 
many as 8 test values (C per card). One 
or more values may be tested for each 
variable indicated on card type D. Addi- 
tional card type E is needed for more 
than 6 values. 

Card Type F 

815 Data are retrieved from this array for as 
many as 8 layer-variable values in addi- 
tion to the ones tested for using card 
types D and E. 

Card Type G 

1615 Data are retrieved from this array for as 

many as 21 drill-hole variables. An ad- 
ditional card, type G, is needed for more 
than 16 test-hole variables. 

Card Type H 

1015 Data are retrieved from this array for as 

many as 10 discrete types of material 
based upon size-distribution analyses. 

Card Type I 

1016 Data are retrieved from this array for as 

many as 10 discrete colors of material. 

Card Type J 

FIO.O Maximum elevation to bo searched (feet 

above mean sea level) when selecting data 
from horizontal subdatum planes. 

PIO.O Minimum elevation to be searched (feet 

above mean sea level) when selecting data 
from horizontal subdatum planes. 

PIO.O The distance in feet to be stepped down for 
subsequent sub datum planes from 

SUBDA to SUBTS. 

Card Type K 

1616 Data are retrieved from this array for as 

many as 16 increments (feet above bed- 
rock) for which data are retrieved. 

Card Type L 

lOFS.O Data are retrieved from this array for as 
many as 10 discrete permeability values. 

Card Type M 

lOFS.O Data are retrieved from this array for as 
many as 10 discrete storage-coefficient 
values. 



Variable 

name 

Option 

Format 

IMO 

t—t 

d 

Card 

16 


IDY 

G,I 

15 

lYR 

G,I 

15 

ITIM 

G.I 

15 

JMO 

G,I 

15 

JDY 

G,I 

15 

JYR 

G.I 

15 

JTIM 

G,I 

15. 


Description of va^^iable 

Type N 

Month of beginning- of time period for which 
ground- water levels are averaged when 
data are retrieved from the saturated 
zone. 

Day of beginning of time period. 

Year of beginning of time period. 

Time (military) of beginning of time 
period. 

Month of ending of time period. 

Day of ending of time period. 

Year of ending of time period. 

Time (military) of ending of time period. 



APPENDIX B.-FORMATS USED FOR DATA STORAGE 


INPUT -CARD FORMAT USED TO STORE DRILL-HOLE DATA 



INPUT -CARD FORMAT USED TO STORE LAYER DATA 



INPUT -CARO FORMAT USED TO STORE WATER-LEVEL DATA 


IRl 

m 

17 



17 

mm 

F7.2 

IT 

16 



* 

d 

z 

,j 

*j 

LJ 

* 

1 CODE NO. 

H ^ S 

S s 

MILITARY TIME 

DEPTH TO WATER 

FEET BELOW 

TOP OF PIPE 


n 

DEPTH 


TIME 

DEPTH 

COLUMN 64-80 
BLANK 


« ALL I FORMAT INPUT DATA IS RIGHT JUSTIFIED 
#* CODE f ■ RECORDER-CHART DATA — CODE 2 •TAPE -DOWN DATA 
NO ZERO PUNCHES NEEDED IN BLANKS ON CARDS 


Figure B-1, — ^Input card formats. 
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APPENDIX C.-SOURCE LISTING FOR STORAGE AND RETRIEVAL 

PROGRAMS 


CARD 

1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

99997 
C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

99998 
C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


APPENDIX C 


RETRIEVAL PROGRAM FOR HYOROGEOLOG I C DATA 
SUBROUTINES USED NAMED (CODE* KTESXi COLOR* SERCH* MATRL* TAPES* JULDY* 
WATLV) 

define file 1 ( 1000*161 *U* ISEC) #2 ( 2 1 • 1000 * U » I SECl ) * 3 ( 1 0 * 3 20 * U * Im ) 
dimension V(21 ) *VV(21) *VARI8) 

DIMENSION MAJTSi 10} ♦MACOLI 10) * WAT EL (1000) » A INC (20) *IPUNH(6) 

COMMON I WELL (1000) *IVAR(0 1 pTVALU(8) *IaOV(8) »WVAR ( 1 000 ) ♦ I STOR ( 2 1 J 
COMMON HEAD (2 >252) * T I TLE ( 2 * 10 8 ) »JWFLL( 1000) *IADVR(8 1 *STRAT(^0*8) 

COMMON ISECf ISECl iISEC2*A<VAL (10) >SVAL(10) * KTEST » <SAT * LTEST 
headers read IN FOR LABELING VARIALBLES* 

CONTROLS FOR SELECTING WELL VARIABLES* 

FORMAT! 1615) 

NWELL-NUMBFR OF WELLS TO BE SEARCHED* 

LVAR -NUMBER OF LAYER VARIABLES TO BE TESTED. 

LVARl-NUMBER OF LAYER TEST VALUES* TVALU (I)> TO BE READ IN. 

VALUES OF (TVALU (I) FORMAT ( 6F12 . 4 t /2F 1 2 * 4 ) ) 

NUM -NUMBER OF ADDITIONAL LAYER VARIABLES NOT USED FOR SEARCHING. 

NVAR -NUMBER OF WELL VARIABLES 


ICONO- 


0 
1 
0 
1 
0 
1 
0 
N 

0 
1 

CONTINUE 


ICOL- 


IDAT- 


INC 


ILOG- 


ILVL- 


IGSUO- 


KTESK- 

KTESS- 


KSAT- 


PUNCH A 


-BYPASS 

-CALLS IN SUBROUTINE FOR TESTING TYPE OF MATERIAL IN THE LAVERS. 
-BY-PASS 

-CALLS IN SUBROUTINE FOR TESTING LAYER COLORS. 

-BY-PASS 

-READS IN SUBSTRATUM ELEVATION TO BE TESTED. (SUBDA) 

-BY PASS 

-READS IN NUMBER OF INCREMENTS IN FEET ABOVE BEDROCK TO BE 
SEARCHED. 

“WRITES OUT LAYER VARIABLE SEARCHED 
-WRITES OUT ALL LAYERS IN A WELL 


BY PASS 

WATER level CALLED FOR SEARCHING IN SATURATED ZONE 
BY PASS 

STEPS DOWN INCREMENTALLY BELOW GROUND SURFACE 
BY PASS 

NUMBER OF PERMEABILITY VALUES TO BE USED 
BY PASS 

NUMBER OF STORAGE COEFFICIENT VALUES TO BE USED 
BY PASS 

K AND S VALUES ARE AVERAGED FOR TOTAL THICKNESS 

IN ILVL AND A IN KSAT TO AVERAGE K AND S VALUES IN THE 

SATURATED ZONE 

FORMAT FOR VARIABLES NWELL THROUGH KSAT IS 1515. 

CONTINUE 

IPUNH (N) - 'O* NO PUNCH U* PUNCH 

THE FOLLOWING PUNCH COMMANDS ARE IN SIX SEPARATE FIVE COLUMN FIELDS. 

THE NUMBER FOLLOWING EACH COMMAND IS THE FIELD NUMBER* 

LAYER DATA (1) 

NET THICKNESS (2) 

WELL VARIABLES (3) 

COLOR OF lithology 14) 

MATERIAL TYPE OF LITHOLOGY (5) 

W6T. AVE. PERM. AND WGT. AVE* STOR. COEFF. (6) 

VALUES OF SUSOA* SUBTS* SUBIN (FORMAT (3F10.0)) 

SUBDA-SUBSTRATUM ELEVATION 
SUBTS-MINIMUM ELEVATION TO BE SEARCHED 
SUBIN-INCREMENT IN FEET TO BE SUBTRACTED FROM SUBDATUM 

SUBROUTINE TAPES LOADS TEST HOLE AND LAYER DATA FROM TAPES TO DISC FILE 
FOR SEARCHING. 

CALL TAPES 
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65 DO 1000 J=lf2 

66 C READ LAYER VARIA8LE MAMES FROM CONTROL CARDS 

67 1000 READ(5»106MHEADIJiI) •1 = 1*1081 

68 106 FORMAT! 7 2A1*/'i36A1) 

69 1M«1 

70 C WRITES LAYER VARIABLE NAMES ON DISC 

71 WRITEO' IMIHEAD 

72 DO 999 J=l*2 

73 C READ TEST HOLE VARIABLE NAMES FROM CONTROL CARDS 

74 999 REA015*107) (HEAD(J*1) *1=1*2521 

75 107 F0RMAT(72A1/72A1/72A1/36A11 

76 1M=3 

77 C WRITES TEST HOLE VARIABLE NAMES ON DISC 

78 WR1TE(3' IMIHEAD 

79 1012 READt5* 100 1 NWELL *LVAR ♦LVARl ♦NUH *NVAR * ICOND* ICOL * IDAT « 1 NC * I LOG * IL 

80 1VL*1GSUB*K.TESK*KTESS*KSAT 

81 C TEST ON NUMBER OF TEST HOLES TO BE SEARCHED 

82 IFINWELLl 1013*1013.4000 

83 C read TEST HOLE NUMBERS TO BE SEARCHED (16151 INTO ARRAY 

84 4000 REAO(5*1001 ( IPUNHI n *1 = 1*61 

85 REA0(5*1001 ( IWELLdl »I = 1*NWELLI 

86 100 FORMAT! 1615) 

87 C TEST ON NUMBER OF LAYER VARIABLES TO BE TESTED 

88 4001 IF(LVARl4003*4003*4002 

89 C READ VARIABLES TO BE TESTED INTO ARRAY 

90 4002 REAO!5* 1001 ( IVAR! I 1 * 1 =1 *LVAR1 

91 C TEST ON NUMBER OF TEST VALUES TO BE TESTED FOR 

92 4003 IFILVAR114005*4005*4004 

93 C READ VALUE OF LAYER VARIABLE TO BE TESTED INTO ARRAY 

94 4004 READ!5*101MTVALU!n *I = 1*LVAR11 

95 101 FORMAT 1 6F 12 *4 */2F 12 .4 1 

96 C TEST ON NUMBER OF ADDITIONAL LAYER VARIABLES NEEDED! 1-71 

97 4005 IF!NUM14007 *4007*4006 

98 C read additional LAYER VARIABLES INTO ARRAY 

99 4006 REA0(5>100) ! lAOVRI 1 1 * 1 = 1 ♦NUM) 

100 C TEST ON NUMBER OF ADDITIONAL TEST HOLE VARIABLES ! 1-21 1 

101 4007 1FINVAR)4009*4009*4008 

102 C READ ADDITIONAL TEST HOLE VARIABLES INTO ARRAY 

103 4008 REA0!5*100) ! ISTOR! I 1 *I=1*NVAR) 

104 C TEST IF SUBDATUM ELEVATION IS TO BE SEARCHED 

105 4009 IF! 10AT)3021*3021*3023 

106 C TEST IF A HEIGHT ABOVE BEDROCK. IS TO BE SEARCHED 

107 3021 IF! INC13022*3022*3024 

108 3024 READ(5*1112)IA1NC(I1 *1»1*INC) 

109 1112 FORMAT! 16F5, 11 

no GO TO 3124 

111 3022 SUSDA>0*0 

112 GO TO 3124 

113 C READ FROM ONE CAR0*!1IELEV TO BE SEARCHED *! 2 ) ELEV TO STOP ON * ! 3 ) 

114 C INCREMENT TO STEP DOWN-FEET 

115 3023 REA015*500ISUBDA*SUBTS*SUelN 

116 500 F0RMAT!3F10,0) 

117 C TEST IF MAJOR OR MINOR LITHOLOGIC TYPE WANTED 

118 3124 IF(1COND)5002*5002*5001 

119 C READ CODE OF LITHOLOGY WANTED INTO ARRAY 

120 5001 REAO!5*100)(MAJTS!1) *l=l*ICOND) 

121 C TEST IF MAJOR OR MINOR COLOR OF LITHOLOGY WANTED 

122 5002 IF! ICOL)5004*5004*5003 

123 5003 READ!5*100) (MACOLd) *I = l*ICOLI 

124 5004 LIST=NUH+LVAR 

125 KTEST»KTESK+KTESS 

126 IF(KTESK)1122*1122*1120 

127 1120 READ(5«112l)(AKVAL!I)*I=l.KTESK) 

128 1121 FORMAT! lOFS.O) 

129 1122 IF(KTESS|1124. 1124, 1123 

130 1123 REAOI5*1121)!SVAL!n *I-1*KTESS) 

131 1124 CALL C0DE!NWELL*LVAR*LVAR1*NUM*NVAR*IC0N0*IC0L.1DAT*INC*IL0G*IPUNH 

132 1*ILVL*MAJTS*MACOL*LIST*SUBOA*SUBTS*SUBIN*ISTOR* lADVR.KSAT* IGSUB.KT 

133 2ESK*KTESS*AKVAL*SVAL*A1NC) 

134 IF(ICONO)5006*5006*5005 

135 5005 CALL MATRL ! I COND .NWELL *MAJTS » I PUNH ! 5 1 1 
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136 5006 IF (I COL I 5080«50a0f 5007 

137 5007 call COLOR ( I COL *rNWFLL *MAC0L , I PUNH ( 4 M 

138 5080 IFI IGSUB|2*2,1 

139 1 SUROA=1.0 

UO AMIN = 9999. 

lAl GO TO 8002 

142 2 IFlSUBOAI8000«8000t8002 

l<t3 8000 IFt ILVL) 8002i8002»800l 

144 8001 CALL WATL V ( NWELL » I WELL * WATEL ) 

145 8002 IFtLIST)6004,6004i30l4 

146 3014 IMnl 

147 C READS HEADERS OF LAYER VARIABLES FROM DISC 

148 RFAOO' IMIHEAO 

149 DO 1003 Kal,2 

150 J=0 

151 C test FOR NUMBER OF VARIABLES TO BE SELECTED 

152 IF(LVAR)903*903t902 

153 902 DO 1002 JC»1,LVAR 

154 C SELECTS HEADERS FOR VARIABLES 

155 ISTRTbI IVARl JC)*12~H) 

156 IENO«lSTRT + n 

157 DO 1001 L“1STRT*IEN0 

168 JKJ+I 

159 1001 TlTLE(K*J)aHEAD(K.LI 

160 1002 CONTINUE 

161 C TEST FOR ADDITIONAL LAYER VARIABLES 

162 IF(NUMH003*1003*903 

163 903 DO 2002 JC«liNUM 

164 C SELECTS HEADERS FOR LAYER VARIABLES 

165 ISTRT“( IA0VR(JC1«12-11) 

166 IEND-ISTRT+11 

167 DO 2001 L»ISTRT»IEND 

168 J^J+l 

169 2001 TITLCtK*J)aHEAO(K>LI 

170 2002 CONTINUE 

171 1003 CONTINUE 

172 C TEST ON SUBDATUM ELEV 

173 IF(SUB0A)5016»5016»3012 

174 5016 IFt ILVL»30l3f3013t5017 

175 C LABELS TEST HOLE N0.» SUBDATUM» DEPTH TO TOP OF LAYER# ELEVATION OF THE 

176 C TOP OF A LAYER* AND ELEVATION OF THE TOP OF SATURATED ZONE* 

177 5017 WRITE(6>109( 

178 WRITB(6»50le) 

179 5018 FORMAT! ' + WN SUB LO L ELEV SATURATED ZONE' I 

180 GO TO 3013 

101 3012 WRITE(6tl09) 

182 C labels TEST HOLE NO.* SUBOATUM* DEPTH TO TOP OF LAYER. ELEV OF TOP OF 

103 C LAYER IN WHICH SUBDATUM IS FOUND 

184 109 FORMAT! 'O' ) 

105 WRITE!6*108) 

186 108 FORMAT! '+'*' WN SUB LD L ELEV' I 

187 3013 IEND=J 

108 DO 1004 K»l*2 

109 C WRITES HEADERS FOR VARIABLES 

190 1004 WRITE!6*115) !TITLEIK*J) *J=1*IENDI 

191 ns FORMAT !' '*22Xt 96A1 ) 

192 ELVl-0.0 

193 ICONT«0 

194 JCONT-0 

195 C SET BRANCH •! 

196 C READ A TEST HOLE NUMBER TO BE SEARCHED 

197 7001 SUB-SUBOA 

190 ICONT-ICONT+1 

199 IF! INC)1110*U10*1U3 

200 1110 IINC*! 

201 GO TO 1114 

202 1113 IlNCsINC 

203 1114 DO nil II«1*IINC 

204 DO 10 1C«1*NWELL 

205 LTEST-O 

206 IBRN»l 


31 



207 

208 

209 

210 
211 
212 
213 
2H 

215 

216 
217 
216 

219 

220 
221 
222 
223 
22 ^ 

225 

226 

227 

228 

229 

230 

231 

232 

233 
23 ^ 

235 

236 

237 

238 

239 

240 

241 

242 

243 


ISEC=21 

JBRN-2 

SUBOA^SUB 

IF(SUBDA)5000i5000f300l 
5000 lF(lLVL)3001*3001i502 
502 JBRN*1 

C READS TOTAL NUMBER OF LAYERS FOR EACH TEST HOLE FROM DISC 

3001 READ (2' ISEOWVAR 
IKKalWELH IC) 

C BRANCHES TO DIFFERENT PARTS OF PROGRAM 
GO T0{3002t3003>53) *IBRN 

C SETS MAXLA EQUAL TO TOTAL NO. OF LAYERS FOR EACH TEST HOLE 

3002 MAXLA=WVAR( IKK) 

GO TO 3004 

C READS GROUND SURFACE ELEV FOR EACH TEST HOLE FROM DISC 

3003 GSELV«WVAR( IKK) 

GSURFaGSELV 

C SETS SECTOR TO FIND BEDROCK ELEV 
ISEC»7 

C SET BRANCH » 3 

IBRN«3 
GO TO 3001 

C READS BEDROCK ELEV FOR EACH TEST HOLE FROM DISC 
53 8EDRK»WVAR( IKK» 

C TEST IF SUBOATUM ABOVE BEDROCK IS DESIRED 
IF( lNCI3006f3006t900l 

C ESTABLISH SUBDATUM ELEV ABOVE BEDROCK 

9001 SU 80 A=>BEDRK+AINCiin 
GO TO 3006 

C TEST IF SUBOATUM ABOVE BEDROCK EXISTS 

3004 IF(SUB0A)9999f9999»3005 
9999 1F( ILVLI9002|9002»3005 

C TEST IF subdatum ABOVE BEDROCK IS DESIRED 

9002 1F( INC)3006i3006i3005 

C SETS SECTOR TO SELECT GROUND SURFACE ELEV# 

3005 ISEC»6 

C SET BRANCH « 2 


244 IBRN*2 

245 GO TO 3001 

246 3006 IFIIGSU8)12001>12001»2000 

247 2000 SUBOAsGSURF 

248 SUB^35UB0A-( IGSUB^^ICONT) 

249 SUBTS*BEDRK 

250 IF (SU8TS» 2003 f 2003 >2005 

251 2003 WRITE(6>2004MWELL( IC) 

252 2004 FORMATS '0* •25( • *NO BEDROCK ELEVATION FOR WELL ' > 1 4 >2 5 ( ’ * M > / / ) 

253 GO TO 10 

254 2005 IF (BEORK-AMIn) 2006# 12001 #12001 

255 2006 AMIN=»8EDRK 

256 12001 ISEC2«IWELLMC) 

257 ILAY«1 

258 READ ( 1 » I 5EC2 ) STRAJ 

259 SUM»0.0 

260 THlK»OiO 

261 7005 1WLL=»IWELL( IC) 

262 C TEST IF TEST HOLE HAS A LOG 

263 7006 1F(MAXLA)10#10#200 

264 C TEST HOLES WITH MORE THAN TWENTY LAYERS ARE WRITTEN ON PRINTER. 

265 200 IFIMAXLA-20)203»203i202 

266 202 WRITE{6#104) IWELLI IC) 

267 104 FORMAT!' '# 'LAYER NUMBER GREATER THAN 20 ON TEST HOLE '#I4) 

268 GO TO 10 

269 C READS THROUGH TEST HOLE LAYERS 

270 203 IF( ILVL)9997f9997f5009 

271 5009 SUBDA»WATEL! IKK) 

272 9997 00 90 JLAY«lfMAXLA 

273 C SUMS LAYER THICKNESS 

274 SUM»SUM + STRATI ILAY#2) 

275 C TEST IF SUBDATUM ELEV IS READ 

276 5 CONTINUE 

277 5008 IF(SUBDA)3008»3008#3007 
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278 3007 IF(SUB0A-GSELVK016. <*016*10 

279 9016 IF(SUBDA-BEDRK) 10*9017*9017 

280 C SUBTRACT LAYER THICK.MESS OFF GROUND SURFACE 

281 9017 ELEV^GSElV-STRAT ( luAY ,2 I 

282 C test new GROUND SURF AGAINST SUBOATUM ElEV 

283 IF ( ELEV-SUBDAI3010*3010*3009 

289 C RESETS GROUND SURFACE TO TOP OF NEXT LAYER 

285 3009 GSELV=ELEV 

286 C INCREASE LAYER NO, BY 1 

287 1LAY=ILAY+1 

288 GO TO 90 

289 C COMPUTES DEPTH TO TOP OF LAYER CONTAINING SUBOATUM 

290 3010 GO TO ( 50 10 * 50 1 2 1 • J8RN 

291 5010 JBRN=2 

292 SUM=GSURF-SUBOA 

293 STRAT(1LAY*2)=SUBDA-ELEV 

299 ELEV=SUBDA 

295 GSELV=SUBOA 

296 SUM=SUM+STRAT< ILAY,2» 

297 IF(KSAT)3007, 3007*7008 

298 7008 CAUL KTESX ( I LAY * S7 , 58 * AKVl * 5V 1 1 

299 5012 DT0LA=GSURF-GSELV 

300 C TEST IF ADDITIONAL LAYER VARIABLES 

301 3020 IF(NUM)3016, 3016*3017 

302 3017 DO 3015 1UJ»1,NUM 

303 C SETS K. TO ADDITIONAL LAYER VARIABLE 

309 K=IAOVR(IJJI 

305 C SETS IV TO THE NUMBER OF THS ADDITIONAL LAYER VARIABLES PLUS THE NUMBER OF 

306 C LAYER VARIABLES TO BE SELECTED 

307 IV-IJU+LVAR 

308 C STORES ADDITIONAL LAYER VARIABLES IN ARRAY 

309 3015 VAR( IV)=STRATnLAY*Xl 

310 IF(LVAR)3016*3025*3016 

311 3016 DO 3011 IJ=1*LVAR 

312 C SETS K TO LAYER VARIABLE TO BE SELECTED 

313 K»IVAR(IJ) 

319 C STORES LAYER VARIABLES SELECTED IN ARRAY 

315 3011 VAR( IJ)>>STRAT( 1LAY*K) 

316 C COMPUTES DEPTH TO TOP OF LAYER WHICH CONTAINS SUBDATUM 

317 3025 ELV»GSURF-(SUM-STRAT(lLAY, 2 n 

318 C SETS IDT TO DEPTH TO LAYER 

319 IDT=DTOLA 

320 C SETS ISUB TO SUBDATUM ELEVATION 

321 ISUB»SUBOA 

322 IF (KSAT ) 300*300* 7003 

323 C TEST IF LOG OF ALL LAYERS IS WANTED 

329 300 IF( ILOG)7050. 7050, 7003 

325 C WRITE DATA OUT FOR ALL LAYERS 

326 7050 1F(1GSUB)7051*7051.7003 

327 C LABELS SUBDATUM ELEVATION* DEPTH TO LAYER* ELEVATION OF TOP OF LAYER 

328 C CONTAINING SUBOATUM* AND SELECTED VARIABLES. 

329 7051 IF(eLV-ELVl)7003*10»7003 

330 7003 WRITE (6*113) IWLL* I SUB * IDT . ELV * ( VAR( 1) . 1 = I * L 1ST ) 

331 113 FORMAT!' ' * 1 3 ♦ 1 5 * 1 9 * F8 . 2 * SF 12 ,9) 

332 C TEST FOR PUNCHED OUTPUT. 

333 IF( IPUNHt 1) )99»99*1050 

339 1050 WRITEI 7* 1051 ) IWLL* ISUB* IDT tELV*t VARt I ) * 1 = 1 (LIST I 

335 1051 FORMAT(I3*15*I9*F8.Z*5F12.9*/»5F12,9) 

336 99 1F(KSATI301*301»93 

337 301 IF( ILVL)10*10»95 

338 95 ILAY=ILAY+1 

339 C TEST LAYER VARIABLE VALUES SEARCHED. 

•>90 C LVARl IS the NUMBER OF VALUES SELECTED FOR EACH LAYER VARIABLE TESTED. 

3A1 3008 IF(LVAR1)9018*9018*909 

392 909 00 6 JJ«=1*LVAR 

393 K=IVAR(JJ) 

399 DO 6 J=l* LVARl 

395 C TEST DATA IN LAYER AGAINST TEST VALUE 

396 IF(TVALU{ J) -STRATI JLAY*K) 16*9011,6 

397 6 CONTINUE 

398 GO TO 9 
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349 c test for additional layer variables 

350 4011 THIK»TH1K+STRAT( JLAY.ai 

351 A018 IF(NUMI6003i6003 »4010 

352 4010 DO 8 IJ=1»NUM 

353 KMAOVRMJI 

354 C SUMS NUMBER OF LAYER VARIABLES TESTED AND ADDITIONAL LAYER VARIABLES 

355 IV=IJ+LVAR 

356 8 VAR( IV) =STRAT( JLAY.K) 

357 C TEST ON NUMBER OF LAYER VARIABLES TO BE SEARCHED 

358 6003 IF(LVAR)6001t6001»6000 

359 6000 DO 7 IJ=1»LVAR 

360 C SETS K TO THE VARIABLE TESTED 

361 K=IVAR(IJ) 

362 C SETS VAR(IJ) TO THE VALUE STORED IN LAYER DATA 

363 7 VAR( I JI»STRAT( JLAY*K) 

364 C TEST ON VARIABLES TO BE SEARCHED FOR OR ADDITIONAL LAYER VARIABLES 

365 6001 IF(LIST)6004»6004*6002 

366 6002 1F( ILVL)6020.6020«9 

367 6020 WRiTE(6*103) IWLL»JLAY» ( VAR( I I *I«!l iLlST I 

368 103 FORMAT!' ' » 1 OX i2 1 5 . 8F12 .4 ) 

369 C TEST FOR PUNCHED OUTPUT 

370 IF! IPUNH! 1 1 )9»9»6006 

371 6006 WRITE17i6007l IWLL*JLAY» (VAR(IJ) »IJ»1*LIST) 

372 6007 F0RMAT!214ieF9.2/»8X»8F9.2/»8X«5F9.2l 

373 C RETURN TO READ NEXT LAYER 

374 9 IF(KTEST)90»90>6008 

375 6008 IF! ILVL)92i92»91 

376 91 KLAY=ILAY-1 

377 GO TO 920 

378 92 XLAY»JLAY 

379 C SUBROUTINE KTESX COMPUTES AVERAGE PERMEABILITY AND STORAGE COEFFICIENT 

380 C VALUES. 

381 920 call KTESX(XLAY*S7iS8fAKVliSVll 

382 90 CONTINUE 

383 93 IF!KTESTI6011i6011.6009 

384 6009 AKV1=AXV1/S7 

385 sviasvi/se 

386 WRITE(6i 6010) IWLLiAKVltSVl 

387 6010 FORMAT!' * • ' HOLE NUMBER* '»I4t/»' WGT. AVE.PERM." '»F10.0>/*' WGTi 

388 lAVE. STORAGE COEFFICIENT* '.F10.3) 

389 IF! IPUNH!6) 16011 »6011»11100 

390 11100 WRITE!7»11101)IWLL»AXVI»SV1 

391 11101 FORMAT! I4.F5,0»F5.3) 

392 C RETURN TO READ NEXT TEST HOLE 

393 6011 IFILVARlll0fl0»4019 

394 4019 WR1TE16 »161 1 IWLL»TH1X» iTVALU! J) tJ*l»LVARl ) 

395 161 FORMAT!' '•' HOLE NO. 'tlSf' NET THICK 'iFlE.O.' CODES '•5X.6F4.0) 

396 IF! IPUNHI2) )10. 10(4020 

397 C PUNCHES DATA ON CARDS 

398 4020 WRITE!7i1611IWLL(THIK((TVALU!J)iJ=1iLVAR1I 

399 10 CONTINUE 

400 1111 CONTINUE 

401 C TEST IF SUBOATUM ELEVATION IS SELECTED 

402 6005 IF!ILVLI99. 99(98 

403 98 SUBDA=0.0 

404 99 IF!SU80A)6004(6004(7000 

405 C TEST IF POINT ABOVE BEDROCK IS SELECTED 

406 7000 IF! INC)9003(9003(6004 

407 C STEP DOWN SUBDATUM TO BE SEARCHED 

408 9003 IF! IGSUB)9004(9004(9005 

409 9004 SUBDA=SUBDA-SUBIN 

410 IF! SUBDA-SUBTS >6004 ( 6004 ( 7002 

411 7002 ELV1*ELV 

412 GO TO 7001 

413 C TEST IF SUBOATUM WANTED IS BELOW LOWER LIMIT 

414 9005 IF!SUBDA“AMIN)6004(6004(7001 

415 C RESET ELFVATION OF TOP OF LAYER 

416 C test if additional TEST HOLE VARIABLES ARE SELECTED 

417 6004 lF!NVAR)1011(10n(ll 

418 C SET DISC SECTOR 

419 11 IM*3 
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4^20 

^21 

422 

423 

424 

C 

425 

426 

C 

427 

C 

428 

1010 

429 

430 

12 

431 

432 

433 

434 

435 

436 

C 

437 

1005 

438 

1006 

439 

440 

441 

1007 

442 

C 

443 

1008 

444 

105 

445 

446 

447 

C 

448 

449 

C 

450 

451 

c 

452 

453 

c 

454 

455 

456 

c 

457 

458 

459 

20 

460 

21 

461 

462 

c 

463 

5400 

464 

465 

C 

466 

1054 

467 

1055 

468 

C 

469 

22 

470 

471 

C 

472 

C 

473 

474 

23 

475 

1009 

476 

477 

1011 

478 

479 

1013 

480 

C 

481 

C 

482 

c 

483 

484 

485 

486 

487 

488 

489 

490 

c 


READ headers FOR ADDITIONAL TEST HOLE VARIABLES 

READ ( 3 ' IMI HEAD 

JSTRTal 

JEND=9 

DO 1009 IL=1»3 

TEST IF NUMBER IS MORE THAN NINE 

IF t NVAR-JENO) 1010, 1010.12 

SETS NUMBER OF VARIABLES TO BE WRITTEN 

JEND=NVAR 

DO 1007 K=l,2 

J = 0 

READS NINE HEADERS FROM DISC STORAGE 

DO 1006 JC= JSTRT »JENO 

ISTRT=( IST0R( JCI#12-1U 

!EN0=ISTRT+11 

DO 1005 L=lSTRT»tEND 

J = J+1 

TITLE(K,JI=HEADIK»L) 

CONTINUE 

CONTINUE 

IENO=J 

DO 1008 K=l,2 

WRITES NINE HEADERS ON EACH PRINTER LINE 
WRITE (6,105) (TITLE ( K » J ) » J= 1 , 1 END ) 

FORMAT!' 'tlOSAl) 

TEST IF SUBOATUM ELEVATION IS SELECTED 

DO 22 J=1»NWELL 

DO 20 I“JSTRT»JEN0 

STORES SELECTED VARIABLES ON DISC 

ISEC1 = IST0R( I I 

READS variables FROM DISC 

READ (2 ' ISECl IWVAR 

SETS IK TO TEST HOLE NUMBER 

IK-slWELLI J) 

SETS IKK TO VARIABLE 
IKKalSTORI I ) 

V( IKK)=WVAR( IK) 

CONTINUE 

DO 21 laJSTRT.JEND 
K»ISTOR ( I ) 

VV( I )=V(K) 

WRITES VARIABLES FOUND AT SUBOATUM ON PRINTER 
WRITE (6,5400) ( VV ( 1) , I = JSTRT , JEND ) 

FORMAT)' ',9F12»3) 

TEST FOR PUNCHED OUTPUT 
IF( IPUNHO) ) 22,22,1054 
WRITe(7,1065) (VV( I ) ,lpJSTRT,JEND) 

FORMAT ( 6F1 2,3, /> 6F12.3 »/ •6F12. 3, /»3F 12.3) 

PUNCHES VARIABLES FOUND AT SUBDATUM, SIX PER CARO 
CONTINUE 

TEST IF ALL variables SELECTED HAVE BEEN WRITTEN 

IF (JEND-NVAR)23, 1011, 1011 

RESET TO NEXT NINE HEADINGS AND VARIABLES 

JSTRT=JEND+1 

JEND=JeN0+9 

CONTINUE 

CONTINUE 

GO TO 1012 

STOP 

END 


subroutine CODE DISPLAYS A LIST OF OPTIONS SELECTED BY THE USER. 
RETRIEVAL CONTROL CARD CODING IS DISPLAYED ON PRINTER FOR EACH JOB 
SUBROUTINE CODE ( NWELL ,LVAR ,LVAR1 ,NUM ,NVAR , I COND , ICOL , lOAT , INC, ILOG 
1 » IPUNH, ILVU,MAJTS»MACOL,1LIST,SU8DA,SUBTS,SUBIN,ISTOR,IAOVK,KSAT,I 
2GSUB,KTESK,KTESS,AKVAL,SVAL,AINC) 

DIMENSION AKVAL( 10) ,SVAL( 10) ,AINC(20) 

DIMENSION MAJTS! 10) ,MACOL ( 10 ) , ISTOR I 2 1 ) , I ADVR (8 ) , I PUNH (6 ) 

COMMON IWELL( 1000) , IVARIB ) ,TVALU(8) 

WRlTe(6,120) 


SETUP. 
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491 120 FORMATCl't' IF TEST CODES ARE EQUAL TO 1. THEN TEST IS USED. IF E 

492 lOUAL TO 0. THEN TEST IS BY-PASSED '♦// I 

493 WRITE(6»100) NWELL 

494 100 FORMAT!' '.'NUMBER OF TEST HOLES TO RE SEARCHED= '.131 

495 WRITE(6.133) ( U'ELLt I 1.1=1 .NWELL) 

496 133 FORMAT!' '.'TEST HOLES THAT ARE BEING TESTED ARE ' . / .4 II 24 1 5 » / I ) 

497 WRITE16.101)LVAR 

498 101 FORMAT! 'O' . 'NUMBER OF LAYER VARIABLES TO BE TESTED= '.12) 

499 IF ILVAR12.2.1 

500 1 WRITEI5.102) ! IVAR! I ) . 1=1 »LVAR) 

501 102 FORMAT!' '.'THE VARIABLES THAT ARE TO BE TESTED ARE'/»8I5) 

502 2 WRITE!6.103)LVAR1 

503 103 FORMAT! 'O' . 'NUMBER OF TEST VALUES* *.I2) 

504 IF!LVARn4.4.3 

505 3 WRITEI6.104)(TVALUII) .I=1.LVAR1) 

506 104 FORMAT!' '.'TEST VALUES OF EACH LAYER ARE ' . / . 8F 10. 4 ) 

507 4 WRITEI6.105)NUM 

508 106 FORMAT! 'O' ♦ 'NUMBER OF ADDITIONAL LAYER VARIABLES* '.12) 

509 IF(NUM)6.6.5 

510 5 WRITEI6.106 ) ! I AOVR! 1 I . 1 = 1 .NUM) 

511 106 FORMAT!' '.'ADDITIONAL LAYER VARIABLES ARE'./iSIS) 

512 6 WRITE!6.116)IL1ST 

513 116 FORMAT! 'O' . 'TOTAL NUMBER OF LAYER VARIABLES THAT ARE CALLED* ’.12) 

514 WR1TE!6.107)NVAR 

515 107 FORMAT! 'O' .'NUMBER OF TEST HOLE VARIABLES* '.12) 

516 1FINVARI8.8.7 

517 7 WRITE!6.108)1ISTORII).I=1.NVAR1 

518 108 FORMAT!' '.'THE TEST HOLE VARIABLES ARE'»/.2U5) 

519 8 WR1TEI6.109) IDAT 

520 109 FORMAT! 'O' . 'SUBDATUM ELEVATION CODE* '.121 

521 IFII0AT)10.10.9 

522 9 WRITEI6.110)SUBDA»SUBTS.SUB1N 

523 110 FORMAT!' '.'ELEVATION TO BC SEARCHED* '.FIO.O./.' ELEVATION TO TER 

524 IMINATE SEARCHING* '.FIO.O./.' STEP-DOWN INCREMENT* '.FIO.O) 

525 10 IF(INC)132.132.130 

526 130 WRITE!6.131)INC.(AINC! I) .I=1.1NCI 

527 131 FORMAT! 'ONUMBER OF HEIGHT INCREMENTS ABOVE BEDROCK TO BE SEARCHED 

528 1= '.12./.' THE INCREMENT VALUES ARE ' . / . 20F 1 0. 2 ) 

529 132 WRITEI6.112)ICONO 

530 112 FORMAT! 'O' . 'TEST IF MAJOR OR MINOR LITHOLOGY TYPE WANTED* '.12) 

531 IF!ICONO)12.12.11 

532 11 WRITE!6.113HMAJTS!I).1»1.IC0NDI 

533 113 FORMAT!' '.'LITHOLOGY TYPE CODES ARE'./. 815) 

634 12 WRITE!6.114) ICOL 

535 114 FORMAT! 'O' . 'TEST IF MAJOR OR MINOR COLOR OF LITHOLOGY WANTED*'. 13) 

536 IF! IC0L114.14.13 

537 13 WRITE!6.115)!MACOL!I). 1=1. ICOL) 

538 115 FORMAT!' '.'LITHOLOGY COLOR COOES ARE'./.aiS) 

539 14 WRITE!6.117) ILVL 

540 117 FORMAT! 'O' . 'TEST IF WATERLCVEL OF SATURATED ZONE IS TO BE USED* '. 

541 1121 

542 WRITE!6.121 IKSAT 

543 121 FORMAT! 'OTEST IF SATURATED ZONE IS TO BE LOGGED* '.121 

544 WRITEI6.1181IL06 

545 118 FORMAT! 'O'. 'TEST IF ALL LAYERS ARE TO BE LOGGED* '.121 

546 IF! IGSUB) 124.124.122 

547 122 WRITEI6.123I1GSUB 

548 123 FORMAT! 'OSTEP DOWN INCREMENT FROM GROUND SURFACE ELEVATION* '.12. 

549 1* FEET’ I 

550 124 IF!KTESK)129.129.125 

551 125 WRITE!6.126)KTESK 

552 126 FORMAT! 'ONUMBER OF K-PERMEABILITY AND S-STORAGE COEF. VALUES* '.12 

553 1) 

554 WR1TE!6.127I!AKVAL!1I .I*1.KTESK) 

555 127 FORMAT!' '.'THE <-PERMEAB I L 1 TY VALUES ARE ' . / . 10F8 . 0 ) 

556 WRITEI6.128I !SVAL! I ) .1=1 .KTESS) 

557 128 FORMAT!' '.'THE S-STORAGE COEF. VALUES ARE ' . / . 1 OF8 .3 I 

558 129 WRITE!6.1191!1PUNH!II .1*1.6) 

559 119 FORMAT! 'O' . 'TEST FOR PUNCHING DATA '»//. 5X .' LAYER OATm* '.I2./.5X.' 

560 INET THICKNESS* ' . 1 2 . / . 5X ' TEST HOLE DATA* ' . 1 2 . / » 6X . ' COLOR OF LITHOLOGY* 

561 1 ' .12./. 5X* 'MATERIAL TYPE OF LITHOLOGY* ’.I2./.5X.' W6T. AVE . PFRM 
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2. AND WGT* AVE« STOR* COEFo 
RETURN 
END 


S12 ) 


562 

563 

564 

565 

566 

567 
563 

569 

570 

571 

572 

573 
5 74 
575 
5 76 

577 

578 

579 
500 
581 

502 

503 

504 

505 

586 

587 
583 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 
601 
602 

603 

604 

605 

606 

607 

608 

609 

610 
611 
612 

613 

614 

615 

616 

617 

618 

619 

620 
621 
622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

632 


C 

c 

C SUBROUTINE KTESX COMPUTES AVERAGE K. AND S VALUES WEIGHTED ON LAYER 
C THICKNESS 

SUBROUTINE KTESXt I LAY » S7 » S8 ♦ AKVl i SV 1 ) 

COMMON IWELLUOOO) >IVAR(8> •TVALU(8) flADVIBI #WVAR(1000) •ISTOR(2l) 
COMMON HEAD (2 >2 52 J » T I T L E I 2 > 108 ) t JWELL* 1000) iEA0VR(8) i STRATI 20 >8) 
COMMON ISEC» ISECl •I5EC2 *AKVAL (101 tSVALl 10) t K TES T »KSA T j lT EST 
lF(LTEST)10*10il5 
10 S7»0*0 
S8«0»0 
AKVl^O.O 
SVlaO*0 
LTEST=^LTEST+l 
15 J7=STRAT( ILAYi7) 

J8«STRAT( ILAY>6 ) 

IF( J7)21 i21»20 

20 S7 = S7>^STRAT( ILAY»2) 

AKV1 = AKV1+(AKVAL( J7)^^STRAT( ILAYi2) ) 

21 IFI J8)26f26i25 

25 S8aS8+STRAT( ILAYf2) 

SVlnSVl+(SVAL( J8)^^5TRAT ( lLAYt2) ) 

26 RETURN 
END 

C 

C 

SUBROUTINE COLOR UCOLR iNWELL •MACOLt 1 PUNH ) 

C SUBROUTINE TO READ AND TEST FOR MAJOR OR MINOR COLOR OF LITHOLOGY 
DIMENSION MACOLI 10) 

COMMON IWELL(IOOO) I IVAR(8) tTVALUIB) » IA0V(8J tWVARtlOOO) i ISTOR(2i) 
COMMON HEAD(2i252) # T I TLE ( 2 • 108 ) i JWELL< 1000) i lAOVRI 8) >STRAT(20*8 J 
COMMON ISECUSECl»rSEC2 

100 FORMAT!' S'HOLE NO. 2X •* LAYER NO. 2X #' MAJOR COLOR ' f 2X ♦ ' M I NOR CO 
ILOR' fAXt 'ELEV ' i3X« 'DEPTH* • 3X ^ ♦ TH I CKNE SS ' fSX.'NET THICK (MAJOR) NE 
2T THICK(MINOR) ' ) 

101 FORMAT I ' ' *2XiI3i8X*I2flOXf I2fr20X#F6.1f2XF6.1f6X#F6.2) 

102 FORMAT! ' ' >2X f 1 3 rQX » I 2 »23X i 1 2 • 7X i F6 . 1 » 2X t F6. 1 i6X >F6* 2 > 

103 FORMAT( ♦+' »0OX»F6.2) 

104 FORMAT! '+' t99X»F6.2) 

C WRITE HEADERS FOR COLORS FOUND IN LITHOLOGY 
WRlTE(6il00) 

C ICOLR IS THE NUMBER OF COLORS SELECTED 
00 6 K«l*ICOLR 

C MACOL(K) IS THE NUMERICAL CODE FOR EACH COLOR SELECTED 
JTESTaMACOL(K) 

C NWELL IS THE NUMBER OF TEST HOLES TO BE SEARCHED 
DO 5 IC^^UNWELL 

C CALLS SUBROUTINE NAMED 'SERCH* 
call SERCH(NLAY#GSELV#tC) 

C SETS SECTOR TO TEST HOLE NO. BEING SEARCHED 
ISEC2=IWELL( IC) 

C READS FROM DISC LAYER DATA 
READ 1 1 ' ISEC2)STRAT 
SUM=0.0 
SUM2=0.0 
SUM3-0.0 
J“1 

C SETS IWLL TO TEST HOLE NO. BEING SEARCHED 

52 IWLL = I WELL ( 1 C 1 

C READS layer DATA FOR ALL LAYERS IN EACH TEST HOLE 

53 DO A JJ»liNLAY 

C SETS TYPE TO COLOR OF LAYER 
50 TYPEsSTRATCJiAI 

C SETS THICK TO THICKNESS OF LAYER BEING SEARCHED 
THICK!»STRAT( J»2 ) 

C SUMS LAYER THICKNESS OF ALL LAYERS AS SEARCHED 
SUM=SUM+STRAT(vl«2t 

C KTEST IS A THREE DIGIT NUMBER (COLOR CODE! OF LAYER 
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633 KTEST* I F I X ( TYPE I / 10 

63A C TEST IF COLOR IN LAYER IS SELECTED COLOR 

635 IF (KTEST-JTEST) 2» 1* 2 

636 C SUMl IS THE THICKNESS TO TOP OF LAYER WHERE COLOR IS FOUND 

637 1 SUM1"SUM-STRAT( Ji2) 

638 SUM2»SUM2+TH1CK 

639 C COMPUTES ELEVATION OF TOP OF LAYER WHERE COLOR IS FOUND 

6A0 ELEVbGSELV-SUMI 

6Al C ICOL IS TWO NUMBER CODE 

6A2 ICOL B type 

6A3 C LABELS HOLE NO. i LAYER N0.» COLOR CODE. ELEV OF TOP OF LAYER. DEPTH TO 
but* C LAYER. LAYER THICKNESS 

646 WRITE(6»10l) IWLL. JJ. ICOL. ELEV.SUMl .THICK 

646 C TEST FOR PUNCHED OUTPUT 

647 1FUPUNHI4.4,10 

648 C PUNCHES HOLE NO.. LAYER NO.. COLOR CODE. ELEV. OF TOP OF LAYER. DEPTH 

649 C TO layer. LAYER THICKNESS ON ONE CARO 

650 10 WRITEn.lllIWLL.JJ. ICOL. ELEV. SUMl. THICK 

651 11 FORMAT!' COLORl ' .3 1 5 .F6 , 1 .2X .F6 . 1 . 6X. F6. 2 ) 

652 GO TO 4 

653 2 KTEST=TYPE”FLOATIKTEST)*10. 

654 C TEST IF COLOR ASKED FOR IS MINOR COLOR 

656 !FIKTEST-JTEST)4.3.4 

656 C CALCULATES DEPTH TO TOP OF LAYER 

657 3 SUM1«SUM-STRAT tJ.2) 

658 SUM3bSUM3+THICK 

659 C CALCULATES ELEV. TOP OF LAYER 

660 ELEV=GSELV-SUM1 

661 C ICOL IS A TWO NUMBER CODE 

662 ICOLbTYPE 

663 C LABELS HOLE NO.. LAYER NO.» COLOR CODE. ELEV OF TOP OF LAYER. DEPTH 

664 C TO LAYER. LAYER THICKNESS 

665 WRITE(6. 102) IWLL. JJ. ICOL. ELEV. SUMl .THICK 

666 IF| IPUNH)4.4.12 

667 C PUNCHES HOLE NO.. LAYER NO.. COLOR CODE. ELEV OF TOP OF LAYER. DEPTH 

668 C TO LAYER. LAYER THICKNESS ON ONE CARD 

669 12 WRITE(7.13)IWLL.JJ. ICOL. ELEV. SUMl. THICK 

670 13 FORMAT! 'COLOR2 ' .3 1 5 .F6 . 1 .2X .F6. 1 .6X . F6. 2 ) 

671 4 J«J+1 

672 IF!SUM2)55. 55.54 

673 54 WRITE!6.103)SUM2 

674 55 IF|SUM3)5.5.56 

675 56 WRITE16.104)SUM3 

676 5 CONTINUE 

677 6 continue 

678 RETURN 

679 END 

680 C 

681 C 

682 SUBROUTINE SERCH INLAY .6SELV . I C ) 

683 C SUBROUTINE READS FROM DISC !1) TOTAL NO. OF LAYERS. 12) GROUND 

664 C SURFACE ELEVATION FOR SELECTED TEST HOLE 

685 COMMON IWELL ! 1000) . IVAR ! 8 ) .TVALU ! 8 I » 1 ADV 1 8 ) . WVAR ! 1000 ) . ISTOR(21 ) 

686 common head 12. 252) . T I TLE ! 2 . 108 ) . JWELL ! 1000 ) .lADVR! 8) .STRATI 20.3 ) 

687 COMMON I SEC . 1 SECl » I SEC2 

688 C SET SECTOR NUMBER TO READ TOTAL NO. OF LAYERS FROM DISC 

689 IBRN>0 

690 ISEC»21 

691 3001 REA012' ISEOWVAR 

692 IKKK*IWELL1 IC) 

693 IF! IBRN) 10.10.3002 

694 10 NLAY»WVAR1 IKKK) 

695 C SET SECTOR NUMBER TO READ GROUND SURFACE ELEV, FROM DISC 

696 ISEC*6 

697 IBRNbI 

698 GO TO 3001 

699 C SETS GSELV TO GROUND SURFACE ELEVATION 

700 3002 GSELVbWVARI IKKK) 

701 RETURN 

702 END 

703 C 
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704 C 

705 SUBROUTINE MATRL UCOND »NWELL iMAJTS t IPUNH I 

706 C SUBROUTINE SELECTS MAJOR AND MINOR MATERIAL TYPE OF 

707 C LITHOLOGY 

708 DIMENSION MAJTS (10) 

709 COMMON IWELL( 1000) .1VAR(8) iTVALU(BI *IADV(B) »WVAR(1000) *IST0RI21) 

710 COMMON H£AO( 2 »252 ) . T I TLE I 2 1 1 08 ) i JWE lL ( 1000 I . IADVR(8) .STRATI 20 t8) 

711 COMMON ISEC.ISEC1.ISEC2 

712 100 FORMAT!' '♦ 'HOLE NO. ' . 2X » ' LA YER NO ' ♦2X »' MAJOR TYPE'.SX. 'MINOR TY 

713 IPE' .SX.'ELEV .AX.'DEPTH' i3X. 'THICKNESS' »5X. 'NET THICK(MAJOR) NET 

714 2THICK(MIN0R) ' ) 

715 101 FORMAT!' ' . 2X * 1 3 *ex * 1 2 • lOX . I 5 * 18X »F6. 1 . 2X >F6. 1 . 6X . F6 • 2 I 

716 102 FORMAT!' ' ♦2X . 1 3 .8X » 12 . 22X » I 5 .6X »F6 • 1 . 2X t F6 . 1 t6X . F6. 2 ) 

717 103 FORMAT! '+' .80X.F6. 2) 

718 104 FORMAT! ' + ' .99X.F6. 2) 

719 WRITEle.lOO) 

720 C ICOND IS THE NUMBER OF MATERIAL TYPES SELECTED. 

721 DO 6 K=l. ICOND 

722 C SETS JTEST EQUAL TO MATERIAL CODE SELECTED (2 DIGIT) 

723 JTEST=>MAJTS!K) 

724 C SEARCH THROUGH EACH SELECTED TEST HOLE 

725 DO 5 IC=1.NWELL 

726 C CALLS SUBROUTINE SERCH 

727 CALL SERCH ! NLAY .GSELV . I C I 

728 ISEC2=IWELL( IC) 

729 C read LAYER DATA FROM DISC FOR EACH HOLE 

730 READ ! 1 ' ISEC2 )STRAT 

731 SUM=0.0 

732 J«1 

733 SUM2=0.0 

734 SUM3“0.0 

735 52 IWLL = 1WELL! IC ) 

736 C SEARCH ALL LAYERS IN EACH HOLE 

737 53 00 4 JJ=1.NLAY 

738 C SET TYPE EQUAL TO MATERIAL CODE FOR EACH LAYER !5 DIGIT) 

739 50 TYPE*STRAT! J.31 

740 C SUM IS THE TOTAL THICKNESS INCLUDING THE LAYER BEING SEARCHED 

741 SUM=SUM+STRAT ! J»2) 

742 C SETS KTEST TO A TWO DIGIT CODE OF MAJOR MATERIAL TYPE 

743 KTEST=1FIX(TYPE)/1000 

744 C THICK IS THE THICKNESS OF LAYER BEING SEARCHED 

745 THICK = STRAT( J.2 ) 

746 C TEST IF MAJOR MATERIAL TYPE IS IN SELECTED LAYER 

747 ID1»TYPE/10000. 

748 IF(ID1-2)1002. 30.1002 

749 30 ITYPE*TYPE/1000. 

750 1D2=1TYPE-! ITYPE/10»10) 

751 IF! 102-2)34.31.34 

752 31 ITEST=102+10 

753 1000 IF!ITEST-JTEST)32ili32 

754 32 ITYPE=TYPE/10« 

755 ID4=ITYPE-!1TYPE/10»10I 

756 1001 lF(ID4-2)35.33.35 

757 33 ITEST=ID4+10 

758 1F(ITEST-JTEST)4.1.4 

759 34 ITYPE=TYPE/lOO. 

760 ID23=ITYPF-( ITYPE/100*100) 

761 IF! ID23-JTEST)1001.1.1001 

762 35 ITYPE»TYPE 

763 ITEST=ITYPC-( ITYPE/100*100) 

764 IF! ITEST-JTEST)4.1,4 

765 1002 IF1KTEST-JTEST)2.1.2 

766 C SUMl IS THE DEPTH TO TOP OF LAYER WHERE MAJOR TYPE IS FOUND 

767 1 SUM1=SUM-STRAT( J.2) 

768 C SUM2 IS THE TOTAL THICKNESS OF LAYERS WITH MAJOR TYPE 

769 SUM2=SUM2+THICK 

770 C ELEV IS ELEVATION OF TOP OF LAYER CONTAINING MAJOR TYPE 

771 eLEV=GSELV-SUMl 

772 C SETS type TO AN INTEGER NUMBER 

773 1TYPE=TYPE 

774 C LABELS HOLE NO.. LAYER NO,. MAJOR MATERIAL TYPE. ELEV. OF TOP OF 
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775 C LAYER. DEPTH TO LAYER. LAYER THICKNESS 

776 WRITE16.101) IWLL* J J . I TYPE « ELEV .SDMI . THI CK 

777 C TEST FOR PUNCHED OUTPUT 

770 IF( 1PUNH1A.4|7 

779 C PUNCHES HOLE NO.. LAYER NO.. MAJOR MATERIAL TYPE. ELEV. OF TOP OF 

7B0 C layer, depth TO LAYER. LAY THICK 

781 7 WRITEI7.8) IWLL. JJ.ITYPE.ELEV.SUMI.THICK 

782 8 FORMAT! 'MAJOR' .316. 2F8.1.F10. 2) 

783 GO TO A 

78A C BREAK DOWN TYPE OF MATERIAL CODE TO TWO DIGITS 

785 2 INUMsTYPE/ 1000. 

786 ITYPE=TYPE 

787 C BREAK DOWN TYPE OF MATERIAL CODE TO MINOR TYPE ONLY 

788 C TEST IF MINOR TYPE OF MATERIAL CODE CALLED IS IN LAYER 

789 INUM=TYPE“IFL0AT I INUM)»1000. I 

790 lAOJ’.INUM/lOO 

791 !F( IADJ-1 121.20.21 

792 20 KTEST=INUM/10 

793 GO TO 22 

79A 21 KTEST=INUM-(1ADJ*100) 

795 22 IF(KTEST-JTEST)A.3.A 

796 C SUMl IS THE DEPTH TO TOP OF LAYER WHERE MINOR TYPE CODE IS FOUND 

797 3 SUM1=SUM-STRAT! J.2I 

798 SUMS.'SUMS+THICK 

799 C ELEV IS THE ELEVATION OF TOP OF LAYER WHERE MINOR TYPE IS FOUND 

800 ELEV»GSELV~SUM1 

801 C LABELS WELL NO.. LAYER NO.. MINOR MATERIAL TYPE. ELEV OF TOP OF 

802 C layer. LAYER THICK 

803 WfiITEt6.102> IWLL. JJ.ITYPE.ELEV.SUMI.THICK 

804 IFI IPUNHI4.A.9 

805 9 WRITE(7.10IIWLL. JJ.ITYPE.ELEV.SUMI.THICK 

806 10 FORMAT! 'MINOR' .316. 2F8.1. Flo. 2) 

807 4 J=J+1 

808 C TEST IF NET THICKNESS OF LAYERS WITH MAJOR MATERIAL IS SELECTED 

809 IF(SUM2 155.55 .54 

810 C WRITE NET THICKNESS OF LAYERS WITH MAJOR MATERIAL TYPE 

811 54 WR1TE(6.103)SUM2 

812 55 IF!SUM3 15 .5.56 

813 56 WRITE! 6 .104)SUM3 

814 5 CONTINUE 

815 6 CONTINUE 

816 RETURN 

817 END 

eie c 

819 C 

820 SUBROUTINE TAPES 

821 C SUBROUTINE TAPES READS TAPES USED IN STORAGE AND RETRIEVAL PROGRAM 

822 C TAPE number 8 FOR WELL VARIABLE DATA. 

823 C TAPE DRIVE NUMBER 9 FOR LAYER DATA, 

824 COMMON I WELL! 1000) . IVAR IS) .TVALU! 8) .lAOV! 8) .WVARI 10001 t I ST0RI21) 

825 COMMON HEAD! 2 .2521 .TITLE I 2 . 108 ) . JWELL ( 1000) . 1 ADVR! 8) .STRAT (20.8) 

826 COMMON ISEC 1 ISEC I . I SEC2 

827 ISECl=l 

828 DO 10 1=1.21 

829 C READS TEST HOLE DATA VARIABLES FROM TAPE 

830 READia.ENDslDWVAR 

831 C WRITES TEST HOLE DATA ON DISC FILE 

832 10 WRITE!2'tSECl)WVAR 

833 11 DC 20 J=1.1000 

834 C READS LAYER DATA VARIABLES FROM TAPE 

835 READ(9»END=301 IM.I .STRAT 

836 C SELECTS LAYER VARIABLES FOR SPECIFIED TEST HOLES 

837 IF! IM)30.30.20 

898 C WRITEILAYER DATA ON DISC) 

839 20 WRITE! I'lMlSTRAT 

840 30 RETURN 

841 END 

842 C 

843 C 

844 C SUBROUTINE JULOY COMPUTES JULIAN DATE FROM CALENDAR DATE AND MILITARY 

845 C TIME. 
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846 subroutine JULDY t IMO * 1 DY * I YR » IHR •! IME » 

847 NYRaIYR/4 

846 NYRsNYR*4 

849 IF (NYR-IYRI 2*lf2 

850 1 IB*=1 

351 GO TO 3 

852 2 IB<»2 

853 3 JDAY=0 

354 DO 10 I»1«IM0 

855 GO TO I 4 » 5 » 6 • 5 »9 f 5 »9 » 5 » 5 » 9 » 5 « 9 ) * I 

856 4 M0«0 

857 GO TO 10 

858 5 MO«31 

859 GO TO 10 

860 6 GO TO (7t8) »IB 

861 7 MO=a9 

862 GO TO 10 

863 8 MO=28 

864 GO TO 10 

865 9 M0=30 

866 10 JOAY=JOAY+MO 

867 J0AY*JDAY+I0Y 

868 DAYaJOAY-1 

869 0AY=-DAY*1440. 

670 HR*( IHR/100I»60 

871 HRlalHR-l ( iHR/lOOUflOO) 

872 TIME=i0AY+HR+HR1 

873 102 RETURN 

874 END 

875 C 

876 C 

877 C SUBROUTINE WATLV COMPUTES GROUND WATER ELEVATION AT SELECTED TEST HOLES 

878 C WHEN DATA FROM THE SATURATED ZONE IS RETRIEVED 

879 SUBROUTINE WATLV (NWELLS * I WELL tWATEL 1 

880 C WATER LEVEL PROGRAM COMPUTES AVERAGE GROuNO WATER ELEVATION FOR A SELECTED 

881 C TIME PERIOD. 

882 DIMENSION IWELL t 1000 ) *WATEL ( 1000 I 

883 REWIND 4 

884 JWELL1=0 

885 C BEGINNING (IMO* IOY» lYR# ITIM) AND ENDING (JMOt JDY» JVRt JTIM) MONTH. 

886 C day. year. MILITARY TIME OF PERIOD FOR WHICH AVERAGE WATER LEVEL IS 

887 C COMPUTED IS READ FROM DATA CARD 

888 1 READ! 5.101 I IMO.IDY. IYR.ITIM.JMO.JDYiJYR.JTIM 

889 101 FORMAT! 1615) 

890 WRITE(6»110) IMO . IDY . I YR . I T IM. JMO . JOY . JYR . JT IM 

891 110 FORMAT! 'O' .25! 'TIME PERIOD FOR SEASONAL ' . 1 2 * ' / ' » 1 2 . ' / ’ . 1 2 . ' / 

892 l'.I4.' To ' . 12. ' . 12 • '/' .12.'/' • I4.1X.25I '«' ) .//) 

893 1IYR»IYR 

894 IF11MO)107.107.5 

895 C SUBROUTINE JULDY COMPUTES TIME IN MINUTES FROM MONTH. DAY. YEAR. ANO 

896 C MILITARY TIME 

897 C BTIM-BEGINNING TIME !JULIAN DATE IN MINUTES) OF PERIOD FOR WHICH WATER 

898 C LEVELS ARE AVERAGED 

899 S CALL JULDY 1 1 MO . I OY . 1 YR . I T IM .BT IM ) 

900 IF! IYR-JYR)2.3.3 

901 C EYEAR-WHEN PERIOD FOR WHICH WATER LEVELS ARE AVERAGED CROSSES END OF YEAR 

902 C THEN EYEAR lEND OF FIRST YEAR. JULIAN DATE IN MINUTES) IS CALUCLATEO 

903 2 CALL JULDY!12»31*IYR.2359»EVEAR) 

904 500 IIYR=IIYR+l 

905 IF! MYR-JYR) 501 .502.502 

906 C STIME-WHEN PERIOD FOR WHICH WATER LEVELS ARE AVERAGED SPANS A FULL YEAR 

907 C THAT YEARS TIME !JULIAN DATE IN MINUTES) IS CALCULATED 

908 501 call JULDYI12.31.IIYR.2359.STIME) 

909 EYCAR=EYEAR+STIME 

910 GO TO 500 

911 C EETIM-THE JULIAN DATE IN MINUTES FROM JANURARY 1 OF THE FINAL YEAR TO THE 

912 C END OF THE SELECTED TIME PERIOD 

913 502 CALL JULDY ! JMO .JOY .JYR .JT IM. EET IM ) 

914 ETlMEsEETIM 

915 eTlM=EETIM+EYEAR 

916 GO TO 6 
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917 C ETIM-ENO OF SELECTED PERIOD (JULIAN DATE IN MINUTES) 

918 3 CALL JULDY( JMO»JOY*jyR»JTIM*ETII^) 

919 ETIME=ETIM 

920 C TTIM-TOTAL TIME (MINUTES) OF SELECTED TIIYE PERIOD 

921 6 TTIM»ETIM“BTIM 

922 C rate of water LEVEL CHANGE IS COMPUTED FROM ELEVATIONS READ FROM TAPE 

9?3 C WITHIN THE SELECTED TIME PERIOD. 

92A DO 1110 I=1iNWELLS 

925 IF( IWFLL( I )-JWELL1)7.1000.8 

926 7 REWIND A 

927 JWELLl'O 

928 8 KKK=IWELL( I ) 

929 JBRN=1 

930 JBN=1 

931 ST1M=0.0 

932 SGH«0.0 

933 TIMZ»0.0 

93A TMoO.O 

935 10 READ(A*END»107) NWELL iNCODE «NYR »NMO »NOY *NT IM t ELEV 

936 IF(NWELL-IWELL(I))10«20.85 

937 20 CALL JULDY ( NMO tMOY.NYR iMT IM iT IME ) 

938 GO T0,(3000»3001) iJBRN 

939 3000 JBRN=2 

9A0 1F( IYR-NYR)86*30*3010 

9A1 3001 lF(lYR-NYR)AOt30»3010 

942 30 IF(DTIM-TIME)3003»3002»3010 

943 3002 WGH»0.0 

944 ELEV2*ELCV 

945 TIMEBsTIME 

946 60 TO 3011 

947 3003 IF(NYR1-NYR)3005 >3004.3004 

948 3004 JBN=2 

949 GO TO 40 

950 3005 JBN=3 

951 GO TO 40 

952 3010 NYR1=NYR 

953 GO TO 10 

954 40 BACKSPACE 4 

955 BACKSPACE 4 

956 READ( 4 >END» 107 )NWeLL>NCODE»NYR>NMO>NOY»NTIM»£LeVl 

957 call JULDYINM0»NDY>NYR>NT1M>T1MA) 

958 45 REAO(4>FNO=107)KWFLL>KCOOE<KYR*KI»»0>KOY»KT1M*ELEV2 

959 IFlKWELL-lWELLd ) )86)46>87 

960 46 1F((KYR-NYR)“1I3014>3013>3012 

961 C A BREAK IN RECORD GREATER THAN ONE YEAR IS PRINTED IN OUTPUT 

962 3012 WRITC(6>109) IWELLI I ) iNMO >NOY >NYR >NT lM>KMO .KDY t KYR > KT IM 

963 109 FORMAT! ' WELL NUMBER '>I4»' BREAK IN RECORD FROM '>3I3>I5t' TO •> 

964 1313(15) 

965 GO TO 1005 

966 3013 IF(KMO+( 12~NMO)-12)30l4>3012*30i2 

967 3014 CALL JULDY ( KMO*KOY >KYR >KT IMi T IMEB ) 

968 GO T0(206(207>3006) (JBN 

969 3006 CALL JULOY ( 12 « 31 >NYR >2359 > TM ) 

970 DIFF«( (TM-T1MA)+BTIM) + (TII^,EB“8T1(^) 

971 DIFF1=TIMEB-BT1M 

972 GO TO 208 

973 206 call JULDY ( 12 >31 >NYR >2359 i TM ) 

974 D1FF=(TM-TIMA)+TIMEB 

975 0IFF1=(TM-6TIM)+TIMEB 

976 GO TO 208 

977 207 DIFF=TIMEB-TIMA 

978 DIFF1=TIME8-BTIM 

979 208 SL0PE=(ELCV2-£L£V1)/0IFF 

980 GHAaELCV2-(SLOPE*DIFFl) 

981 0IFF=DIFF1 

982 STIM=STIMtDlFF 

983 IFISTIM“TTIM)1002>1011tl010 

984 1010 01FF=D1FF-(TM-ETIM+TIMEB) 

935 GHB=ELEV2--SL0PE*(TM-ETIM+TIMEB) 

986 GO TO 1012 

987 1011 'GHB=ELEV2-SL0PE*(TIMEB-ET-1M) 
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988 1012 WGH=( ( (GHA+GHB)-»t5MDIFF)/TTIM 

989 GO TO 9001 

990 1002 &HB=>EL£V2 

991 WQH=( ( (GHA+GHe)*.5)*0IFF)/TT[M 

992 3011 SGH=SGH+WGH 

993 GHA=ELEV2 

99A T1MZ=TIMEB 

995 50 REA0( A»END»107)NWELLtNCODE*NYRtNMO*NOYtNTIM*QHB 

996 IF(NWELL-IWELl.< I ) )B6»51 *87 

997 51 CALL JULDY I NMO *N0Y * NYR *NT IM* T IMX ) 

998 IF(TIMZ“TIMX)201*201*200 

999 200 CALL JULDY ( 12 * 3 1 iNYR-1 * 2359 * TM ) 

1000 DI FF= ( ( TM-TIMZ )+TlMX) 

1001 GO TO 203 

1002 201 DIFF^^TIMX-TIMZ 

1003 203 STIM=STIM+DIFF 

1004 IF(STIM“TT IMI56i800*70 

1005 55 WGH= ( (( (GHB+GHAI*.5)*D1FF (/TTIM) 

1006 SGH=SGH+WGH 

1007 GHA=GHB 

1008 TIM2=TIMX 

1009 GO TO 50 

1010 800 WGH=( ( ( (GH8+GHAI*.5»*0IFF)/TT1MI 

1011 GHA=6HB 

1012 TIMZ-TIHX 

1013 GO TO 9001 

1014 70 BACKSPACE 4 

1015 80 BACKSPACE 4 

1016 READ(4*ENO=107 )NWELL,NCODE*NYR.NMO»NOY.NTtM*EL£Vl 

1017 IF(NWELL-IWELL(n)86*81*87 

1018 81 GHA=ELEV1 

1019 CALL JULOY(NMO*NDY.NYR*NTIMtTIMEA) 

1020 READ(4*EN0b 1O7 )NWELL*NCODe*NYRl*NMO»NOY*NTIM*ELEV2 

1021 CALL JULDY(NM0*NDY»NYR1 *NTIM*TIMEB) 

1022 IF(NYRl-MYRie05*805*801 

1023 801 CALL JULDY ( 1 2 * 3 1 »NYR *2359 tTMl ) 

1024 DIFF=(TH1-TIMEA)+TIMEB 

1025 SL0PF*(ELEV2-ELEV1)/DIFF 

1026 IF(JYR-NYR)e03*803*802 

1027 802 GHB«ELEV2-(SL0PE*(TmEB-ETlME) > 

1028 WGHb ( ( ( (GHB+GHA)*. 5 l»t ( TMl-TIMEA) +ET 1 ME I ) /TTIM ) 

1029 GO TO 9001 

1030 803 GHB«ELEVl+iSLOPE»<ETlME-TIMEA) > 

1031 GO TO 9000 

1032 805 SL0PE=( ELEV2-ELEVUZtTIMEB-TIMEA) 

1033 6HBoELEV2- (SLOPE* (TIMEB-E TIME) ) 

1034 IF (ETIME-TIMEA) 9000 *9009* 9000 

1035 9009 ETIME=TIME0 

1036 9000 WGHn( ( ( (GH9+GHA)**5)*(ET1ME-TIMEA) )/TTIMI 

1037 9001 WGHP=SGH+WGH 

1038 GO TO 1006 

1039 C WHEN NO RECORD IS FOUND FOR SELECTED TIME PERIOD THIS IS PRINTED IN OUTPUT 

1040 85 WRITE(6 *103) IWELL( t ) * IMO* IDY* lYR* ITIMi JMO* JDY*JYR* JTIM 

1041 103 FORMAT!' '*' WELL NO# '*14*' NO RECORD FOR PERIOD Of '*313*15*' TO 

1042 1 ’*313*15) 

1043 GO TO 1005 

1044 C when selected TIME PERIOD STARTS BEFORE RECORDS ON TAPE THIS IS PRINTED IN 

1046 C OUTPUT* 

1046 86 WRITEI6*104| IWELL(I) 

1047 104 FORMAT!' '*' WELL MO* '*14*' PERIOD OF RECORD STARTS BEFORE ACTUAL 

1048 IRECORO') 

1049 GO TO 1005 

1050 C WHEN selected TIME PERIOD GOES BEYOND RECORDS ON TAPE THIS IS PRINTED IN 

1051 C OUTPUT 

1052 87 WRITE!6.105) IWELLH ) 

1053 105 FORMAT!' '*' WELL NO* '*14*' PERIOD OF RECORD BEYOND ACTUAL RECORD 

1054 I'l 

1055 C NUMBER 99999 STORED IN WATEL (KKK) INDICATES A MISSING RECORD FOR THESE 

1056 C SELECTED TIME PERIODS 

1057 1005 WGHP=99999. 

1058 1006 WATEL! KKK )»WGHP 
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1059 1000 JWELL1=IWELI.( n 

1060 1110 COMUNUE 

1061 107 RETURN 

1062 END 

1063 C 
lOe^t C 

1065 C END OF RETRIEVAL PROGRAwSi STORAGE PROGRAMS FOLLOW. 

1066 C 

1067 C «*♦* STORAGE PROGRAM FOR TEST HOLE DATA 

1066 DEFINE FRE 1 1 2 1 > 1000 *U . IM ) 

1069 DIMENSION WELLl I 1000 ) .WELL I A ) *NN ( 4) 

1070 DIMENSION ISTART(21) 

1071 DO 105 1=1*21 

1072 105 ISTART<I)=0 

1073 C READ FROM CARD THE NUMBER OF TEST HOLE VARIABLES TO BE STORED ON DISC 

1074 fiEAO(5«10O) MOVAR 

1075 100 FORMAT! 141 

1076 C READ A VARIABLE NUMBER* FOUR VARIABLE VALUES* AND FOUR TEST HOLE 

1077 C NUMBERS FROM CARDS 

1078 99 READ(5*200)NVAR * INNI I ) *1 = 1*4) *(WELL( 1) * 1 = 1 *4) 

1079 200 FORMAT! I2*4I5*4FIU0) 

1080 IF!. WAR) 4 *4,1 

1081 1 IM^NVAR 

1082 IFIISTARTINVARI.EO.O) GO TO 6 

1063 REAOII'IM) WELLl 

1084 JH=IM-1 

1085 GO TO 8 

1086 6 DO 7 1=1*1000 

1087 7 WELLl! I 1«0. 

1088 8 00 3 t*l»4 

1089 IFINN! I ) )3 *3*2 

1090 2 IW=NNin 

1091 WELLll IW)=WELL( r ) 

1092 3 CONTINUE 

1093 C WRITE TEST HOLE VARIABLE ON DISC 

1094 WRITE! I'lMl WELLl 

1095 rSTART!NVAR) =1 

1096 GO TO 99 

1097 C WRITE TEST HOLE VARIABLE DATA FROM DISC TO TAPE 

1098 4 IM»1 

1099 00 5 J=l*NOVAR 

1100 READd'IHI WELLl 

1101 WRITE! 9) WELLl 

U02 C WRITE variable NUMBER* TEST HOLE NUMBER* ANO VARIABLE VALUE'ON PRINTER 

1103 WR1TE!6*101 » I * ( J iWELLl 1 J ) * J = 1 * 1 000 1 

1104 101 FORMAT! 'O' • 'variable NUMBER * • 1 2 ♦ / / * 7 ! 1 4 * 2X * FI 1 . 0 ) ) 

1105 5 CONTI, NUE 

1106 WRITFI6*102) 

1107 102 FORMAT t' I'* 'TEST HOLE VARIABLE DATA HAS BEEN STORED ON TAPE') 

1108 STOP 

1109 END 

1110 C 

1111 C, 

1112 C 

1113 C *»** STORAGE PROGRAM FOR LAVER DATA **** 

1114 define file U 1000,161 *U*IMJ 

1115 DIMENSION STRAT(20,8I *NWELLS(1000) 

1116 C NOWEL IS THE NUMBER OF TEST HOLES FOR WHICH DATA IS STORED 

1117 R£AD!5ilOOINOWEL 

1118 100 FORMAT! I4) 

1119 ICOUN=0 

1120 DO 3 1C«1*N0WEL 

1121 C NWELU IS THE TEST HOLE NUMBER 

1122 C STRAT !l.JI IS THE LAVER VARIABLES STORED IN EACH LAYER 

1123 C NOLAY IS THE NUMBER OF THE LAYER FOR WHICH VARIABLES ARE STORED. 

1124 READ! 5 » 101 INWELL*! STRAT! l.J) ,J«l*a) *N0LAY 

1125 101 FORMAT! I4,2F4.0*1X,F5 .0,2F3.0*2F7.0*F10.6* 14) 

1126 IF1NWELLI200*200*300 

1127 300 ICOUN=ICOUN+l 

1128 IF(flOLAY-ll3*2»l 

1129 1 READ! 5* 106)1 (STRAT! I *J) *J»l,e) *I*2*NOLAY) 
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1110 106 FORMAT('tX»2FA.OtlX>F5.0*2F3.0»2F7.0iF10.6) 

1131 2 WRITEIl'NWELCINOLAY.STRAT 

1132 NWELLSI IC)=NWECL 

1133 3 CONTINUE 

113A C SEQUENCE THE WELLS IN ASCENDING ORDER 

1135 200 NUM=ICOUN-l 

1136 DO 5 J=ltICOUN 

1137 DO 5 I=1«NUM 

1138 IF(NWELLS( I ) -NWELLS ( I+l ) I5»5f4 

1139 4 HOLD=NWELLSt I ) 

1140 NWELLS(I)»NWELLS(I+1) 

1141 NWELLSl I+1)=HOLO 

1142 5 CONTINUE 

1143 C WRITF LAYER DATA ON TAPE AND PRINTER 

1144 DO 7 IC=l»ICOUN 

1145 IM-'NWELLSI ICI 

1146 READ! 1 ' IMlNOLAYfSTRAT 

1147 IU=IM-1 

1148 WRITE(8) lUtNOLAYtSTRAT 

1149 DO 6 I=l»NOLAY 

1151 103 FORMAT!' ' 1 2 1 5 »8F10. 1 ) 

1152 7 CONTINUE 

1153 IMaO 

1154 WRlTEial IM»NOLAY»STRAT 

1155 WRITE(6»104I 

1156 104 FORMAT! ' 1 ' t 'LAYER DATA STORED ON TAPE'! 

1157 STOP 

1158 END 

1159 C 

1160 C 

1161 C **** WATER LEVEL DATA STORAGE PROGRAM *»»* 

1162 C CALLS SUBROUTINE ITESTEL) 

1163 C PROGRAM STORES PUNCH CARO TAPE DOWN DATA AS WATER LEVEL ELEVATIONS ON 

1164 C MAGNETIC TAPE* 

1165 C master tape OF WATER LEVEL ELEVATIONS IS CONSTRUCTED USING PIPE ELEVATIONS 

1166 C AND tape DOWN DATA STORED ON CARDS OR TAPES. 

1167 C WATER LEVEL DATA FOR SELECTED TEST HOLES CAN BE UPDATED USING THIS 

1168 C PROGRAM* 

1169 C TAPE PARAMETERS READ IN BY DATA STATEMENT 

1170 C WORDS-NUM8ER OF WORDS WRITTEN PER RECORD 

1171 C LRECL-LENGTH OF RECORD SIZE 

1172 C BLKSI2-0LOCK SIZE PARAMETER 

1173 C BPI-eiTS PER INCH FOR OUTPUT TAPE 

1174 C RCAP-RECORD GAPE BETWEEN RECORDS I INCHES) 

1175 REAL44 LRECHINPBK 

1176 DIMENSION KWL 1 4 ) .ELVX 1 4 ) . ELVN 1 1000 1 » I TEST ! 1000 ) » JWELL ! 3 ) 

1177 DIMENSION IWN! 100) , IMOl 100) . lOYI 100 1 » I YR! 100) tELV! 100 1 

1178 DIMENSION I TAPE I 3 ) .NMO ( 3 ) .NDY ( 3 ) »NYR 1 3 ) .NT IM! 3 ) »EL ! 3 ) 

1179 COMMON NYRtNMO»NDY.NTIM.ELt lYR.lMO* IDYfELV.ELVN.ITEST.KTAPE. 

1180 *NCODE. ICOUNT»IPA6E. IWN 

1181 DATA WORDS»LRECL«BLKSIZ.BP1*RGAP/8..24.»2400.»1600.».76/' 

1182 DATA JWELL/340/ 

1183 DO 1 1=1.1000 

1184 ITESTII)=1 

1185 1 ELVN! I )=0.0 

1106 C READ IN TOP OF PIPE ELEVATIONS FOR ALL TEST HOLES. 

1187 C IV-LAST CARD TEST 

1188 C XWL ID-TEST HOLE NUMBERS* 

1189 C ELVX !l)-TOP OF PIPE ELEVATIONS AT EACH TEST HOLE. 

1190 5 READ!5.1000)IVi(KWLlI)»I«1.4l.(ELVX(l)»i=1.41 

1191 1000 FORMAT! 12. 4I5.4F11. 2) 

1192 IFIIV.LE.O) GO TO 15 

1193 00 10 1=1.4 

1194 IFIKWLI I I *LE.O) GO TO 10 

1196 KKK=KWL!1) 

1196 ELVN!K.KK)=.ELVX( I) 

1197 10 CONTINUE 

1198 GO TO 6 

1199 C READ IN PIPE ELEVATIONS CHANGES 

1200 15 1=1 

1201 C IWN ID-TEST HOLE NUMBER 
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1202 

1203 

1204 

1205 
1205 
1207 

120a 

1209 

1210 
1211 
1212 

1213 

1214 

1215 

1216 
1217 
1216 

1219 

1220 
1221 
1222 

1223 

1224 

1225 

1226 

1227 

1228 

1229 

1230 

1231 

1232 

1233 

1234 

1235 

1236 

1237 

1238 

1239 

1240 

1241 

1242 

1243 

1244 

1245 

1246 

1247 

1248 

1249 

1250 

1251 

1252 

1253 

1254 

1255 

1256 

1257 
1256 

1259 

1260 
1261 
1262 

1263 

1264 

1265 

1266 

1267 

1268 

1269 

1270 

1271 

1272 


lYR in IS THE MONTH* DAY* AND YEAR OF PIPE ELEVATION 


C IMO n 1 . IDY (1 I * 

C (ELV ID) CHANGES. 

16 READ! 5* 1001 ) IWN( I > » IMOl I I * lOY I 11 * lYR t 1) *EUVI I ) 

1001 F0RMAT<4I5*F10.2 ) 

IF( IWN< I ) .LE.OI GO TO 17 
ITESTI IWNII I )=0 
1 = 1 + 1 
GO TO 16 

17 M0BS1=I 

C NTAPE-NUMBER OF TAPE DRIVE UNITS USED. 

C ITAPE <n-TAPE DRIVE UNIT NUMBERS. 

C <TAPE-CREATEO OR UPDATED TAPE UNIT NUMBER. 

REA0(5. 1002 IMTAPE*! I TAPE! I ) . I =lfMTApE) .KTAPE 

1002 FORMAT (515 I 

C READ WELL NUMBER FROM TAPES FOR I N I T I L 1 2AT I ON . 

1005 00 1006 IT=1 ,mTAPE 
ITP=ITAPE( 1T> 

R£A0nTP(lO3> JWELLI IT) 

1006 CONTINUE 

1007 WRITE I 6 *9999) I JWELLI I ) * I = 1*NTAPE ) 

.9999 FORMAT I '0*«*»* '*3110*' ##*#+*/) 

C FOLLOV/ING TEST STATEMENTS SELECT DATA FROM TWO OR THREE DATA SETS AND 

C ARRANGES THEM IN ORDER OF ASCENDING TEST HOLE NUMBERS FOR CREATING OR 
C UPDATING DATA TAPE. 

IFIJWFLLIU .EQ 
IFlJWFLLin.EO 
IF( JWFLLd ) .EO 
1F( JWELLI 1 ) .GT 
IF( JWELLl 1 ) .EG 
IFIJWFLLID.EO 
IF( JWELLI 1 ) .EQ 
IFIJWCLLI2) .EQ 
IFIJVELLI2) .EO 
IF(JV)ELLI2) .EQ 
ir< JWELLI3 > .EQ 
IFIJWELLI3I .EO 
tF( JWELLI3 I .EQ 
IF( JV)ELLI 1 I .EO 
IF( JWELLI 1 1 .EO 
tF( JWELLI 1 I .LT 
IF ( JWELLI 1) ,LT 
IFt JWELLI 1 1 .GT 
IF ( JWELLI D.GT 
IFUWELLID.LT 
IFIJWFLLIl) .LT 
IF (JWELLI 1) .GT 
IFIJWELLID.GT 
WRlTE{6tl009) 

1008 FORMATI 'O' * 'THERE 
GO TO 2000 

1009 DO 34 IT=1*NTAPE 

iTP«iTAPEH n 
IFIITAPEIITI.EQ.OI 
BACKSPACE ITP 
NWFLL1=JKELLHTP) 

20 IFIlTP.EQ.l) GO TO 


.AND.JWELL 

(2 

) .EQ.O.AND. Jh'FLL 1 3) • 

EO.OIGO 

TO 

2000 

.and.jwfll 

(2 

) .EQ.O.AND. JWELLI 3) . 

GT.OIGO 

To 

1060 

.AND. JWFLL 

12 

) .GT.O. AND.JWFLL! 3) . 

EO.OIGO 

TO 

1050 

.AND. JWELLI 

(2 

) .EQ.O.AND. JWELLI 3) . 

EQ.OIGO 

TO 

1030 

.AND. JWELL i 

(2 

) .EQ. JWELLI 3) ) 

GO 

TO 

1040 



.AND.JWELL! 

(2 

) .LT. JWELLI 3) ) 

GO 

TO 

1050 



•AND.JWELL! 

(2 

) .GT. JWELLI 3) ) 

GO 

TO 

1060 



•AND. JWELLI 

!1 

) .EQ. JWELLI 3) ) 

GO 

TO 

1020 



.AND. JWELLI 

[1 

I.LT. JWEULI3I ) 

GO 

TO 

1030 



.and. JWELLI 

:i 

) .CiT. JWELLI 3) ) 

GO 

TO 

1060 



.AND.JWELL! 

:i 

).EO.JWELL(2) ) 

GO 

TO 

1010 



.AND.JWFLL! 

:i 

) .LT. JWELLI 2 1 ) 

GO 

TO 

1030 



.AND.JWELL! 

1 

) .GT.JWELLia) ) 

GO 

TO 

1050 




JWELLI 2 I .AND. JV/ELL I 2 ) .EQ. JWELU ( 3 ) 1 GO TO 1009 
JWELLI 2 I .AND.JWELL (2 ) .LT. JWFLL<3) I GO TO 1010 
JWELLI 2 ) .AND. JWELL 1 1 ) »EQ. JWELH3 I ) GO TO 1020 
JWELL I 2 ) .AND. JWELL I 2 ) .EQ. JWELLI 3 ) ) GO TO 1030 
JWELL I 2 ) .AND. JWELL ( 2 ) .EQ. JWELL ( 3 ) ) GO TO 1040 
JWELL (2) .AND. JWELLI 2 ) .LT. JWELL (3 ) ) 60 TO 1050 
JWELLI 2 ) .AND. JWELH2) .LT. JWELL 13 ) ) GO TO 1030 
JWELLI2) .AND.JWELL(2) .GT.JWELLI3> ) GO TO 1030 
JWELLI 2) .AND.JWELL (2 ) .LT. JWELLI 3 ) I GO TO 1050 
JWELLI2I .AN0.JWELLI21 .GT.JWELLO) ) GO TO 1060 

ARE NO TEST COMBINATIONS FOR THESE TEST HOLES.'/) 


GO TO 34 


1 


103 

21 

104 

22 

25 

30 


21 

READ(|TP,103,END=30)NWELL,NCODE. (NMOII ) .NOYII ) *NYRI I ) ,NTIMI I ) .ELI 1 
IREAD=3 

FORMAT! I4*12.3(1X*312*15*F7.2)I 
GO TO 22 

FORMATn4l4I2n4°F8!2')'^^^’^^°°^*^’^°'^* > . NT IM 1 1 ) .EL I 1 ) 

IREAD=1 

IFINWFLLI «N£ .NWELL ) QO TO 25 
CALL TESTELINWELL.NOBSI * IREAD) 

NWELLl=NW£LL 
GO TO 20 

JWELLI ITP)=NWELL 
GO TO 34 
JWELU n TP 1=0 
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1273 ITAPEtm=0 

127A 34 CONTINUE 

1275 GO TO 1007 

1276 1010 DO 45 IT=1»2 

1277 ITP=1TAPE( IT) 

1278 IFUTAPE( IT).EQ»0) GO TO 45 

1279 BACKSPACE ITP 

1280 NWEUL1 = JWELL( ITP) 

1261 35 IF(tTP.EO.l) GO TO 36 

12 82 READ! 1 TP ♦ 103 » EN0=41 »NWELL *NCOOE* ( NMO ( I ) •NDY( I )tNy?!(I)tNTIMtl)*ELtI 

1283 11.1=1.3) 

1284 IREAD=3 

1285 GO TO 37 

1286 36 READ! ITP. 104 . END=4 1 (NWEL L . NCODE* NMO( 1 1 .NDY< 1) .NYR(l).NTIM(l)»EL(ll 

1287 IREAD=1 

1288 37 IF(NWELLI.NE.NWELL) GO TO 40 

1289 CALL TESTELINWELL.NOBSl.lREAOl 

1290 NWELL1=NWELL 

1291 GO TO 35 

1292 40 JWELLI ITP)=NWtLL 

1293 GO TO 45 

1294 41 JWEUL( ITP)=0 

1295 rTAPE(lT)=0 

1296 45 CONTINUE 

1297 GO TO 1007 

1298 1020 DO 65 17=1.2 

1299 GO TO (50.51) .IT 

1300 50 ITP=ITAPE(1) 

1301 IF( ITAPEI 1) .EQ.O)GO TO 65 

1302 GO TO 52 

1303 51 ITP=ITAPE(3) 

1304 IF(ITAPE(3).EQ,0)GO TO 65 

1305 52 BACKSPACE ITP 

1306 NWELLl=JweLL( ITP) 

1307 55 IF(ITP.EO.l) GO TO 56 

1308 READ( ITP. 103.END»61)NW£LL.NC0DE» (NM0( I ) .NOYd ) .NYR( I ) .NT IM( 11 lEL ( I 

1309 11.1=1.3) 

1310 IREAO=3 

1311 GO TO 57 

1312 56 READ! ITP. 104. END»6l )NWELL.NCOOEfNMO( 1 ) »NDY( 1) ».NYR( 1 l.NTIMtD.ELd) 

1313 IREAO=l 

1314 57 IF(NWFLLl.NE.NWELL) GO TO 60 

1315 CALL TESTELINWELL.NOBSl.lREAOl 

1316 NWELL1=NWELL 

1317 GO TO 55 

1318 60 JWELLI ITP)=NWELL 

1319 GO TO 65 

1320 61 JWELLI ITP)=0 

1321 ITAPE(1TP)»0 

1322 65 CONTINUE 

1323 60 TO 1007 

1324 1030 ITP»ITAPE(1) 

1325 IF! ITAPEI 1 ) .EQ.O) GO TO 1007 

1326 BACKSPACE ITP 

1327 NWELH = JWELL(1TP) 

1328 70 IFIlTP.EOtll GO TO 71 

13 29 READ! I TP . 103 .END»80 INWELL .NCODE. ( NMO I J 1 »NOY 1 1 1 .NYRC I 1 » NT IMI 1 1 . EL M 

1330 1).I=1.3I 

1331 IREA0=3 

1332 GO TO 72 

1333 71 READ! ITP. 104. END = 80)NWELL. NCODE .NMOIl) .NDY< 1) >NYRI 1) . NT I M 1 1 ) . EL 1 1 ) 

1334 IREAD=1 

1335 72 IFINWFLLl.NE.NWELLI GO TO 75 

1336 CALL TESTELINWELL.NOBSI.IREAD) 

1337 NWELL1=NWELL 

1338 GO TO 70 

1339 75 JWELL(ITP»=NWELL 

1340 GO TO 1007 

1341 80 JWELLI 1TP)=0 

1342 1TAPE(1I=0 

1343 CO TO 1007 
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1344 

1345 

1346 

1347 

1348 

1349 
13&0 

1351 

1352 

1353 

1354 

1355 

1356 

1357 

1358 

1359 

1360 

1361 

1362 

1363 

1364 

1365 

1366 

1367 

1368 

1369 

1370 

1371 

1372 

1373 

1374 

1375 

1376 

1377 

1378 

1379 

1380 

1381 

1382 

1383 

1384 

1385 

1386 

1387 

1388 

1389 

1390 

1391 

1392 

1393 

1394 

1395 

1396 

1397 

1398 

1399 

1400 

1401 

1402 

1403 


1040 DO 100 IT*^lf2 

GO TO(85»86)iIT 

85 ITP=ITAPE{2) 

IF( ITAPE(2) «£0*0»GO TO 100 
GO TO 87 

86 ITP»ITAPE<3) 

IFnTAPE(3)«EQiO)GO TO 100 

87 BACKSPACE ITP 
NWELLl»JWELUnTP) 

90 IFdTPiEQ.l) GO TO 91 

READ ( ITP f 103 »EN0 = 96)NWELLf NCODEt (NMO( I ) iNOY^ I ) tNYR( I ) »NT 1 M U ) f EL U 
1) 

1READ»3 
CO TO 92 

91 READ< ITPf 104 tEN0=96 INWELL fNCODE^NMOt 1) »NDY( 1 1 iNYRl 1) • N T 1 M ( 1 ) » EL t 1 I 
IREA0=1 

92 IF{NKELLltNE*NW£LLI GO TO 95 
call TESTEL(NWELLiN0BS1 I IREAD) 

NV*ELL1«NWELL 

GO TO 90 

95 JWELL( ITP)=NWELL 
GO TO 100 

96 JWELU(ITP)«0 
ITAPE( ITP)«0 

100 CONTINUE 
GO TO 1007 
1050 1TP»ITAPE(2) 

IF(ITAPE(2)»EOtO)GO TO 115 
BACKSPACE ITP 
NWELLl^JWELLIITP) 

105 IFilTPiEOill GO TO 106 

READnTP»103*EN0«linNWELLiNC0DEf INMOI 1 ) fNDYC 1 1 iNYRI I ) >NT IM ( I ) >EL ( 

in »i«it3} 

IREADa3 
GO TO 107 

106 REAO( ITPi 104iEND»lll INWELLfNCOOE •NMOl 1 ) iNDY ( 1 ) fNYR ( 1 ) iNT IM( 1 I t EL ( 1 
*1 

IREADal 

107 IFINWELLltN£*NWELL) GO TO 110 
CALL TESTELINWELLtNOBSlflREAD) 

NWELL1=NWELL 

GO TO 105 

110 JWELL(ITP}=NWELL 
GO TO 115 

111 JWELL(ITPI*0 
ITAPE( ITP)«0 

115 CONTINUE 
GO TO 1007 
1060 ITP=ITAPEI3) 

IF(ITAPEI3)fEQiO)GO TO 130 
BACKSPACE ITP 
NWELL1»JWELL(ITP) 

120 IF(ITPiEOtl) GO TO 121 

READ I ITP • 103 iEND*l26 )NWFLL f NCODE » (NMO( I ) tNOY tI)fiNYR(I)iNTlM{nfEL( 
II 1 ii^ns) 

IREA0='3 
GO TO 122 

12l^READnTPtX04iEN0«l26)NWELLfNC0DE>NM0( 1) fNOYU ) >NYR( 1 ) tNT I M ( 1 ) t EL 1 1 
IREAD*^! 


1404 

1405 

1406 

1407 

1408 

1409 

1410 

1411 

1412 

1413 

1414 


122 IF(NWELUUNEiNWELL) GO TO 125 
CALL TESTELlNWELLtNOBSl • IREAD) 

NWELLl^NWELL 
GO TO 120 

125 JWELLIITPI^NWELL 
GO TO 130 

126 JWELL(ITPI=0 
ITAPEi ITP)=»0 

130 CONTINUE 
GO TO 1007 

C COMPUTES THE NUMBER OF FEET OF TAPE USED TO STORE WATER LEVEL ELEVATIONS 
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1^^15 C AND THE NUMBER OF RECORDS STORED. 

1416 2000 RECD=ICOUNT 

1417 RECDS=BLKSI2/LRECL 

1416 INPBK=( (W0RDS4BLKSIZ)/BPI )+RGAP 

1419 FTs( <R£CD/RECDS)<tlNP3K)/12. 

1420 WRITE(6*2001I I COUNT .KT APE • FT 

1421 2001 FORMAT( 'O' « NUMBER OF OBSERVATIONS WRITTEN ON TAPE '.I2» 

1422 ft' 5$$$$' /»6X. 'FEET USED TO STORE DATA ON THIS TAPE IS '.F6.1I 

1423 REWIND KTAP£ 

1424 STOP 

1425 END 

1426 C 

1427 C 

1428 C SUBROUTINE TESTEL CALCULATES WATER LEVEL ELEVATIONS AND WRITES ON UPDATED 

1429 C TAPE. 

1430 SUBROUTINE TESTEL I NWELL . NOBSl • I READ I 

1431 COMMON NYRf 3) »NM0(3) iNDY (31 .NTIM ( 3 ) » EL I 3 I » I YR ( 1 00 ) » 1 MO 1 1 00 > > I DY ( 10 

1432 10) lELVI 100) .ELVNdOOO) 1 1 TEST) 1000 I .K.TAPE .NCODE . ICOUNT . I PAGE ♦ I WN ( 10 

1433 10) 

1434 DO 28 1=1 » IREAD 

1435 IF( IREAD. eO.l ) GO TO 1000 

1436 IFINYRI I) .EQ.O) GO TO 28 

1437 IFIITeST(NWELL).EQ.l) GO TO 27 

1438 DO 25 0=1 iNOBSl 

1439 C TEST FOR PIPE ELEVATION CHANGE 

1440 IFINYRI I ) .GE. lYRI Jl.ANO.NMOI 1 1 .GE . IMO( Jl .ANO.NDYd ) .GE. IDYI J) .AND. 

1441 INWELL.EQ. IWN( J) ) GO TO 26 

1442 25 CONTINUE 

1443 GO TO 27 

1444 C NEW PIPE ELEVATION 

1445 26 ELVNINWELL) »ELV( J) 

1446 1TEST(NWELLI=1 

1447 C WATER LEVEL ELEVATION CALCULATED. 

1448 27 FLI I )=eLVN(NWELL)-EL( I ) 

1449 C OUTPUT WRITTEN ON UPDATED TAPE 

1450 1000 V/RtTE(KTAPe*104)NWELL.NCOOE*NM0( I ) »NOYl I ) .NYRI I ) >NT I M( I ) .EL ( I) 

1451 104 FORMAT! I4*4I2i IA.F8. 2 ) 

1452 28 CONTINUE 

1453 RETURN 

1454 ENO 

1455 C END APPENDIX C 
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APPENDIX D.-SOURCE LISTING FOR HYDRAULIC-COEFFICIENT 

RANGE-ASSIGNMENT PROGRAM 


1 C APPENDIX D 

2 C 

3 C 

A C HYDRAULIC COEFFICIENT RANGE ASSIGNMENT PROGRAM 

5 C »»*THIS PROGRAM ASSIGNS HYDRAULIC COEFFICIENT RANGE TO EACH LAYER*** 

6 C THIS PROGRAM REQUIRES A SCRATCH TAPE FOR ASSIGNING CODE VALUES TO LAYER 

7 C DATA 

e DIMENSION Al(35) ♦A 2 (l 5 ) »A3( 10) *AA115I *STRAT(20*e) iA5( 15) 

9 C **THE FOLLOWING CARO READS IN THE NUMBER OF MATERIAL TYPES IN EACH RANGE** 

10 READ15*200) INO»JNO,KNOfLNO*MNO 

11 200 FORMAT(5I5l 

12 C »»*THE FOLLOWING CARD PUTS RANGE 1 MATERIAL INTO AN ARRAY FOR LOADING*** 

13 READ(5«101)(Al(I)»I°ltINO) 

lA C •**THE following CARD PUTS RANGE 2 MATERIAL INTO AN ARRAY FOR LOADING#** 

15 REAOISilOll IA2(I)*I*l»JNO) 

16 C *»»THE FOLLOWING CARO PUTS RANGE 3 MATERIAL INTO AN ARRAY FOR LOADING*** 

17 READ(5»101)(A3(I).I=liKNO) 

18 C »**THE FOLLOWING CARD PUTS RANGE 4 MATERIAL INTO AN ARRAY FOR LOADING*** 

19 REAOI5«101)(A4(1).I=1»LNO) 

20 C ***THE FOLLOWING CARD PUTS RANGE 5 MATERIAL INTO AN ARRAY FOR LOADING*** 

21 RFAO(5*101 1 (A51 1 I »I»l»MNO) 

22 101 FORMAT! I0F8.0) 

23 C **»TH£ FOLLOWING CARO READS THE TEST HOLE NO.» NO. OF LAYERS AND MATERIAL 


26 

27 

2a 

29 

30 

31 

32 

33 
3k 
35 
3 6 
3 / 

38 

39 

40 

41 

42 

43 


^ 10 mu ^END=100UNWELLiNOLAY,STRAT 

DO 1000 I=l#20 
IF (NOUAYtEQiOJGO TO 1002 

C SERIES OF TESTS MATCHES MATERIAL TYPE TO THOSE IN EACH RANGE*^^i^ 

C ***WHEN MATERIAL TYPES MATCH THE RANGE IS ASSIGNED AND STORED ON TAPE^^^^^t 
00 20 J=l»INO 

1FISTRAT(I»3)-A1IJ)) 20»60»20 
20 CONTINUE 

DO 30 JaliJNO 

IFISTRATI I #3I-A2I J) 1 30>70f30 
30 CONTINUE 

00 40 J^liKNO 

IFISTRATI I i3)-A3i J) ) 40#e0i40 
40 CONTINUE 

DO 50 J«liLN0 

IF I STRATI I fSI-AAlJ) ) 50 #90 *50 
50 CONTINUE 

DO 500 Js^liMNO 

IF (STRATI I >3 J-A5IJH500*501*500 


44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 


500 CONTINUE 
GO TO 1000 

C CODE l«Oi 2i0# 3i0* 4i0* OR 5*0 IS ASSIGNED TO 
C RANGE l.Oi 2.0* 3.0* 4.0* OR 5.0 RESPECTIVELY 
C range CRITERIA IS DETAILED IN TEXT 
60 CODE=L*0 
GO TO 100 
70 C0DE«2.0 
GO TO 100 
80 COOE:=3.0 
GO TO 100 
90 CODE’»4.0 
GO TO 100 

501 CODE*»5.0 

100 00 110 K=6>8 
no STRATI I jlO^'COOE 

1000 CONTINUE 

1002 WRITE<9)NWELL*N0LAY*STRAT 
GO TO 10 

1001 REWIND 8 
ENOFILE 9 
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65 rewind 9 

66 100^ READ ( 9 »EN0=1003 I NWELL »NOLAY iSTRAT 

67 WRITE(8)NWELL»NOLAY*STRAT 

68 GO TO 1006 

69 1003 REWIND 8 

70 WRITe(6»10051 

71 c ***the following statement indicates ranges have been stored on tape*** 

72 1005 FORMAT (' ODATA ON TAPE') 

73 STOP 

76 END 

75 C END APPENDIX D 


appendix e.-flow charts and list of variables for 

PREPARING GRAPHICAL PRESENTATIONS 


Flow Charts 


MAP VIEW 


CROSS SECTION 




SCALE RARAMEURS 
XLENG.VLENO.GTEST.BTFSr 
{SEE VARIABLE UST) 



retrieve layer variables 

AND STORE ON DISK , NWELL, 

THICK. KVALU 

{SEE VARIABLE LIST) 

J 

1 

RETRIEVE TEST-HOLE 
VAWABUS a STORE ON DISK 
NWaL.GSURF, TDEPH. MAXU 
(SEE VARIABLE LIST) 


f 

CALL SCALE PLOT 
SUBROUTINE 



CALL GRID AXIS SUBROyTI^E 
AND ANNOTATE THE GRID 


r 

PLOT INDIVIDUAL TEST-HOLE 
LOO AND WATER SURFACE 


Figure E-1.— Flow charts for preparing map views and cross sections. 


51 














List of Variables 


Variable 

name 

Format 

Desoiptwn of variable 


Boundary Card for Cross Section 

XLENG 

PIO.O 

Distance in inches between 

YLENG 

PIO.O 

test-hole-log plots. 

Length in inches of the cross- 

GTEST 

PIO.O 

section plot. 

Reference ground-surface 

BTEST 

PIO.O 

datum from which logs are 
measured on the plot. 
Lower reference bedrock 


datum beneath which test- 
hole logs do not extend, 


Lithology Card for Cross Section 


NWELL 

14 

Test-hole numbers of logs 
used in cross section, 

IVAR 

14 

Layer number for lithology 
of test holes in cross sec- 
tions. 

THICK 

F9.1 

Thickness of lithologic layers 
in feet. 

KVALU 

F9.1 

Hydraulic-coefficient range 
(1-6) of the lithology in 
each layer. 


Test-Hole Card for Cross Section 

NWELL 

18 

Test-hole number for each 
log in cross section. 

GSURF 

P12.3 

Ground-surface elevation 
(feet above mean sea 
level) for the top of each 
log. 

TDEPH 

F12.3 

Total depth of each test hole. 

MAXLA 

12 

Number of lithology layers in 
each test hole in cross sec- 
tion. 

Groimd-Water-Elevation Card for Cross Section 

HWELL 

16 

Test-hole number. 

HCODE 

12 

Code indicating the method 
of measurement of gi^oimd- 
water level (1 — continuous 
recorder, 2— manual tape 
down). 

NYR 

13 

Year when water surface was 
recorded. 

NMO 

13 

Month when water surface 
was recorded. 

NDY 

13 

Day when water surface was 
recorded. 

NTIM 

16 

Time (military) when water 
surface was recorded. 

ELEV 

F8.2 

Ground-water-surface eleva- 


tion (feet above mean sea 
level) in the test hole. One 
card per test hole. 


Variable „ . . . . , 

format Description of vaiiable 

Location Card for Map Views 

Note: four maps are plotted: distributive map, sub- 
datum map, hydraulic-coefficient map, isopachous 
map. 

NW F12.3 Test-hole number for each 

hole plotted in map view. 

ALON F12.3 Longitude in degrees, min- 

utes, seconds of test-hole 
locations. 

ALAT F12.3 Latitude in degrees, minutes, 

seconds of test-hole loca- 

tions. 

Lithology Card for Subdatum Map 

NW 13 Test-hole number of each 

test hole that penetrates 
the specific datum plane. 

ISUB 16 Elevation of the datum plane 

to be plotted (feet above 
mean sea level). 

IDE PH 14 Depth below ground surface 

to the top of the lithologic 
layer in which specified 
datum plane lies. 

ELTOP F8.2 Elevation of the top of the 

lithology layer that con- 
tains the specified datum 
plane (feet above mean sea 
level) . 

RANGE F12.4 Range of hydraulic coeffi- 

cients (1-5) at the speci- 
fied datum. 

THICK F12.4 Thickness in feet of the lith- 

ology layer found at the 
specified datum. 

DAPLT F12,4 This variable name is as- 

signed to the range of K, 
St or T to be plotted. 

Test-Hole Card for Map View 

Note; four maps are plotted: distributive map, sub- 
datum map, hydraulic-coefficient map, isopachous 
map, 

NW 16 Test-hole number. 

THVAR P6.0 Variable THVAR may be 

any of the 21 test-hole 
variables used for plotting. 

Net-Thickness Card for Isopachous Map 

WELL NO, AlO Descriptive literal field. 

NW 16 Test-hole number. 

NET THICK AlO Descriptive literal field. 

THICK F13.0 Net thickness in feet of 

lithologic layers having the 
selected hydraulic-coeffi- 
cient ranges. 
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Variable 

name 


Foyynaf Description of vai iahle 


Variable 

name 


Format Desoiption of variable 


Net- thickness Card for Isopachous Map— Continued 
CODES A6 Descriptive literal field. 

CODE (I) 6F4,0 Nimiber (1-5) of the rang'es 

of hydraulic coefficients, 
that IS, Ky S, or to be 
plotted. 

Card for Ilydrauhc-Coefficient Map 
NW 14 Test-hole number, 


Card for Hydraulic-Coefficient Map — Continued 

AKVAL F6.0 Permeability averaged for all 

layers in the test hole, gpd/ 
ft- of aquifer, The size of 
the fields containing hy- 
draulic coefficients may be 
altcied to accommodate re- 
trieved data. 

ASVAL F6.0 Storage coefficient, dimen- 

sionless value. 

ATVAL Fll.O Transmissibility, gpd /ft of 

aquifer. 
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mmm card for cross sectioh test-hqie card for cross section 


APPENDIX F.-INPUT CARDS FOR GRAPHICAL PRESENTATIONS 
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Figure F-l. — ^Input cards for plotting cross sections. 
















LOCATION CARD FDR HAP VIEW 



PiGUtiB P-2 —Input cards for plotting map views. 






















FifiURE F-3, — Input-card setup for plotting cross sections. 



Figure F-4. — In put- card setup for plotting test-hole data (see figure 10) . 



Figure F-6. — Input-card setup for plotting data at a specified subdatum 

(see figure 12) . 



Figure P-6. — Input-card setup for plotting average peimeability coefficients 

(see figure 14), 



Figure F-*?. — Input card setup for isopachous maps (see figure 16) . 
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